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|. QCD AND NUCLEAR FORCES 


Within the Standard Model of particle physics, the 
strong interactions are described by Quantum Chromo- 
dynmics (QCD). QCD is a fascinating theory with many 
intriguing manifestations. Its structure and interactions 
are governed by a local non-abelian gauge symmetry, 
namely SU(3)color- Its fundamental degrees of freedom, 
the quarks (the matter fields) and gluons (the force carri- 
ers), have never been observed in isolation (confinement). 
The strong coupling constant ag exhibits a very pro- 
nounced running and is of order one in the typical energy 
scales of nuclear physics. The bound states made from 
the basic constituents are the hadrons, the strongly in- 
teracting particles. The particle spectrum shows certain 
regularities that can be traced back to the flavor sym- 
metries related to the fermions building up these states. 
More precisely, there are six quark flavors. These can be 
grouped into two very different sectors. While the light 
quarks (u, d,s) are almost massless and thus have to be 
treated relativistically, bound states made from heavy 
quarks allow for a precise non-relativistic treatment. In 
what follows, we will only consider the light quarks at 
low energies, where perturbation theory in ag is inappli- 
cable (this regime is frequently called “strong QCD”). A 
further manifestation of strong QCD is the appearance of 
nuclei, shallow bound states composed of protons, neu- 
trons, pions or strange particles like hyperons. The re- 
sulting nuclear forces that are responsible for the nuclear 
binding are residual color forces, much like the van der 
Waals forces between neutral molecules. It is the aim of 
this article to provide the link between QCD and its sym- 
metries, in particular the spontaneously and explicitely 


broken chiral symmetry, and the nuclear forces which will 
allow to put nuclear physics on firm theoretical grounds 
and also gives rise to a very accurate calculational scheme 
for nuclear forces and the properties of nuclei. 

This review is organized as follows: In this section, we 
briefly discuss some of the concepts underlying the chiral 
effective field theory of the nuclear forces and make con- 
tact to ab initio lattice simulations of two-baryon systems 
as well as to more phenomenological approaches. Sec. IT 
deals with the foundations and applications of nuclear 
EFT and should be considered the central piece of this 
review. In particular, tests of these forces in few-nucleon 
systems are discussed. Attempts to tackle nuclear matter 
and finite nuclei are considered in sec. III. We end with 
a short summary and outlook. 


A. Chiral symmetry 


First, we must discuss chiral symmetry in the context 
of QCD. Chromodynamics is a non-abelian SU(3)color 
gauge theory with Ny = 6 flavors of quarks, three of them 
being light (u,d,s) and the other three heavy (c, ),t). 
Here, light and heavy refers to a typical hadronic scale of 
about 1 GeV. In what follows, we consider light quarks 
only (the heavy quarks are to be considered as decou- 
pled). The QCD Lagrangian reads 


1 ane s 
Lecp = roe (GuvG’’) + Giy*Dig—GMq 
= Bop - Ma, (1.1) 


where we have absorbed the gauge coupling in the defi- 
nition of the gluon field and color indices are suppressed. 
The three-component vector q collects the quark fields, 
q' (x) = (u(s), d(x), s(x)). As far as the strong inter- 
actions are concerned, the different quarks u,d,s have 
identical properties, except for their masses. The quark 
masses are free parameters in QCD - the theory can be 
formulated for any value of the quark masses. In fact, 
light quark QCD can be well approximated by a fictitious 
world of massless quarks, denoted Loop in Eq. (1.1). Re- 
markably, this theory contains no adjustable parameter - 
the gauge coupling g merely sets the scale for the renor- 
malization group invariant scale Aqgcp. Furthermore, in 
the massless world left- and right-handed quarks are com- 
pletely decoupled. The Lagrangian of massless QCD is 
invariant under separate unitary global transformations 
of the left- and right-hand quark fields, the so-called chi- 
ral rotations, qr > Vrqr ,Vr € U(3) ,f = L, R, leading 
to 3? = 9 conserved left- and 9 conserved right-handed 
currents by virtue of Noether’s theorem. These can be 
expressed in terms of vector (V = L+R) and axial-vector 
(A = L — R) currents 


: ost 
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Ab = Gy"95¢, Ab=Gy"r=¢@, (1.2) 


Here, a = 1,...8, and the Aq are Gell-Mann’s SU(3) fla- 
vor matrices. The singlet axial current is anomalous, and 
thus not conserved. The actual symmetry group of mass- 
less QCD is generated by the charges of the conserved 
currents, it is Gp = SU(3)r x SU(3), x U(A)v. The 
U(1)y subgroup of Go generates conserved baryon num- 
ber since the isosinglet vector current counts the number 
of quarks minus antiquarks in a hadron. The remain- 
ing group SU(3)r x SU(3)z is often referred to as chiral 
SU(3). Note that one also considers the light wu and d 
quarks only (with the strange quark mass fixed at its 
physical value), in that case, one speaks of chiral SU (2) 
and must replace the generators in Eq. (1.2) by the Pauli- 
matrices. Let us mention that QCD is also invariant 
under the discrete symmetries of parity (P), charge con- 
jugation (C) and time reversal (JT). Although interesting 
in itself, we do not consider strong C’P violation and the 
related 0-term in what follows, see e.g. (1). 


The chiral symmetry is a symmetry of the Lagrangian 
of QCD but not of the ground state or the particle 
spectrum — to describe the strong interactions in na- 
ture, it is crucial that chiral symmetry is spontaneously 
broken. This can be most easily seen from the fact 
that hadrons do not appear in parity doublets. If chi- 
ral symmetry were exact, from any hadron one could 
generate by virtue of an axial transformation another 
state of exactly the same quantum numbers except of 
opposite parity. The spontaneous symmetry breaking 
leads to the formation of a quark condensate in the 
vacuum (0|gq|0) = (O|¢@r¢r + Grqz|0), thus connect- 
ing the left- with the right-handed quarks. In the ab- 
sence of quark masses this expectation value is flavor- 
independent: (O|tu|0) = (0|dd|0) = (0|gq|0). More pre- 
cisely, the vacuum is only invariant under the subgroup 
of vector rotations times the baryon number current, 
Hy = SU(3)y x U(1)y. This is the generally accepted 
picture that is supported by general arguments (2) as 
well as lattice simulations of QCD (for a recent study, 
see (3)). In fact, the vacuum expectation value of the 
quark condensate is only one of the many possible order 
parameters characterizing the spontaneous symmetry vi- 
olation - all operators that share the invariance properties 
of the vacuum qualify as order parameters. The quark 
condensate nevertheless enjoys a special role, it can be 
shown to be related to the density of small eigenvalues 
of the QCD Dirac operator (see (4) and more recent dis- 
cussions in (5; 6)), lima—o(0|Gq|0) = —7 p(0). For free 
fields, p(A) ~ \3 near A = 0. Only if the eigenvalues 
accumulate near zero, one obtains a non-vanishing con- 
densate. This scenario is indeed supported by lattice 
simulations and many model studies involving topologi- 
cal objects like instantons or monopoles. 


Before discussing the implications of spontaneous sym- 
metry breaking for QCD, we briefly remind the reader 
of Goldstone’s theorem (7; 8): to every generator of a 
spontaneously broken symmetry corresponds a massless 
excitation of the vacuum. These states are the Gold- 
stone bosons, collectively denoted as pions 7(x) in what 


follows. Through the corresponding symmetry current 
the Goldstone bosons couple directly to the vacuum, 

(0|A°(0)|z) #0. (1.3) 
In fact, the non-vanishing of this matrix element is a nec- 
essary and sufficient condition for spontaneous symme- 
try breaking. In QCD, we have eight (three) Goldstone 
bosons for SU(3) (SU(2)) with spin zero and negative 
parity — the latter property is a consequence that these 
Goldstone bosons are generated by applying the axial 
charges on the vacuum. The dimensionful scale associ- 
ated with the matrix element Eq. (1.3) is the pion decay 
constant (in the chiral limit) 

(0| Af, (0)|n°(p)) = 10 F py , (1.4) 
which is a fundamental mass scale of low-energy QCD. In 
the world of massless quarks, the value of F' differs from 
the physical value by terms proportional to the quark 
masses, to be introduced later, F, = F[1+ O(M)]. The 
physical value of F, is 92.4 MeV, determined from pion 
decay, 7 > VL. 

Of course, in QCD the quark masses are not exactly 
zero. The quark mass term leads to the so-called explicit 
chiral symmetry breaking. Consequently, the vector and 
axial-vector currents are no longer conserved (with the 
exception of the baryon number current) 


OVE = FGM dala, nA = 511M, da} o50- 
(1.5) 
However, the consequences of the spontaneous symme- 
try violation can still be analyzed systematically because 
the quark masses are small. QCD possesses what is called 
an approximate chiral symmetry. In that case, the mass 
spectrum of the unperturbed Hamiltonian and the one 
including the quark masses can not be significantly dif- 
ferent. Stated differently, the effects of the explicit sym- 
metry breaking can be analyzed in perturbation theory. 
As a consequence, QCD has a remarkable mass gap - the 
pions (and, to a lesser extent, the kaons and the eta) are 
much lighter than all other hadrons. To be more specific, 
consider chiral SU(2). The second formula of Eq. (1.5) 
is nothing but a Ward-identity (WI) that relates the ax- 
ial current A“ = dy"y5u with the pseudoscalar density 
P = diysu, 
Op AY = (Mu+ma)P. (1.6) 
Taking on-shell pion matrix elements of this WI, one ar- 
rives at 
M2 = (mu + ma) 


TT 


(1.7) 


where the coupling G, is given by (0|P(0)|(p)) = G,. 
This equation leads to some intriguing consequences: In 
the chiral limit, the pion mass is exactly zero - in ac- 
cordance with Goldstone’s theorem. More precisely, the 


ratio G,/F; is a constant in the chiral limit and the pion 
mass grows as ,/mz + mq if the quark masses are turned 
on. 

There is even further symmetry related to the quark 
mass term. It is observed that hadrons appear in isospin 
multiplets, characterized by very tiny splittings of the or- 
der of a few MeV. These are generated by the small quark 
mass difference m, — mq and also by electromagnetic ef- 
fects of the same size (with the notable exception of the 
charged to neutral pion mass difference that is almost en- 
tirely of electromagnetic origin). This can be made more 
precise: For m, = ma, QCD is invariant under SU (2) 
isospin transformations: q — q’ = Uq, with U a unitary 
matrix. In this limit, up and down quarks can not be dis- 
entangled as far as the strong interactions are concerned. 
Rewriting of the QCD quark mass term allows to make 
the strong isospin violation explicit: 


Heep = m,tut+madd (1.8) 
~ mel (au + da) + a au — dd) , 


where the first (second) term is an isoscalar (isovector). 
Extending these considerations to SU(3), one arrives at 
the eighfold way of Gell-Mann and Ne’eman that played a 
decisive role in our understanding of the quark structure 
of the hadrons. The SU(3) flavor symmetry is also an 
approximate one, but the breaking is much stronger than 
it is the case for isospin. From this, one can directly infer 
that the quark mass difference m, — mq must be much 
bigger than mq — My. 

The consequences of these broken symmetries can be 
analyzed systematically in a suitably tailored effective 
field theory (EFT), as discussed in more detail below. 
At this point, it is important to stress that the chiral 
symmetry of QCD plays a crucial role in determining the 
longest ranged parts of the nuclear force, which, as we 
will show, is given by Goldstone boson exchange between 
two and more nucleons. This was already stressed long 
ago, see e.g (9) (and references therein) but only with 
the powerful machinery of chiral effective field theory this 
connection could be worked out model-independently, as 
we will show in what follows. 


B. Scales in nuclear physics 


To appreciate the complexity related to a theoretical 
description of the nuclear forces, it is most instructive to 
briefly discuss the pertinent scales arising in this prob- 
lem. This can most easily be visualized by looking at 
the phenomenological central potential between two nu- 
cleons, as it appears e.g. in meson-theoretical approaches 
to the nuclear force, see Fig. 1. The longest range part 
of the interaction is the one-pion exchange (OPE) that 
is firmly rooted in QCD’s chiral symmetry. Thus, the 
corresponding natural scale of the nuclear force problem 
is the Compton wavelength of the pion 


Ne = 1/M, ~ 1.5 fm, (1.9) 
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FIG. 1 Schematic plot of the central nucleon-nucleon poten- 
tial. The longest range contribution is the one-pion-exchange, 
the intermediate range attraction is described by two-pion 
exchanges and other shorter ranged contributions. At even 
shorter distances, the NN interaction is strongly repulsive. 


where M, = 139.57 MeV is the charged pion mass. The 
central intermediate range attraction is given by 27 ex- 
change (and shorter ranged physics). Finally, the wave- 
functions of two nucleons do not like to overlap, which 
is reflected in a short-range repulsion that can e.g. be 
modelled by vector meson exchange. From such consid- 
erations, one would naively expect to be able to describe 
nuclear binding in terms of energy scales of the order of 
the pion mass. However, the true binding energies of the 
nuclei are given by much smaller energy scales, between 
1 to 9 MeV per nucleon. Another measure for the shal- 
low nuclear binding is the so called binding-momentum 
y. In the deuteron, y = /~mBp ~ 45MeV « M,, with 
m = 938.2 MeV the nucleon mass and Bp = 2.224 MeV 
the deuteron binding energy. The small value of y signals 
the appearance of energy/momentum scales much below 
the pion mass. The most dramatic reflection of the com- 
plexity of the nuclear force problem are the values of the 
S-wave neutron-proton scattering lengths, 


|a(So)| = 23.8fm >> 1/M,, a(?S,) =5.4fm > 1/M,. 

(1.10) 
Thus, to properly set up an effective field theory for the 
forces between two (or more) nucleons, it is mandatory to 
deal with these very different energy scales. If one were to 
treat the large S-wave scattering lengths perturbatively, 
the range of the corresponding EFT would be restricted 
to momenta below pmax ~ 1/|a(+So)| ~ 8 MeV. To over- 
come this barrier, one must generate the small binding 
energy scales by a non-perturbative resummation. This 
can e.g. be done in a theory without explicit pion de- 
grees of freedom, the so-called pion-less EFT. In such an 
approach, the limiting hard scale is the pion mass. To 
go further, one must include the pions explicitely, as it is 
done in the pion-full or chiral nuclear EFT. The relation 
between these different approaches is schematically dis- 
played in Fig. 2. A different and more formal argument 
that shows the breakdown of a perturbative treatment 
of the EFT with two or more nucleons is related to the 
pinch singularities in the two-pion exchange diagram in 
the static limit as will be discussed later in the context 
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FIG. 2 Scales in the two-nucleon problem and the range of 
validity of the corresponding EFTs as explained in the text. 
Here A, is the hard scale related to spontaneous chiral sym- 
metry breaking, with A, ~ M,, with M, = 770 MeV the mass 
of the rho meson. 


of the explicit construction of the chiral nuclear EFT. 

In addition, if one extends the considerations to heav- 
ier nuclei or even nuclear matter, the many-body sys- 
tem exhibits yet another scale, the Fermi momentum kr, 
with kp ~ 2M, at nuclear matter saturation density. 
This new scale must be included in a properly modified 
EFT for the nuclear many-body problem which is not 
a straightforward exercise as we will show below. It is 
therefore not astonishing that the theory for heavier nu- 
clei is still in a much less developed stage that the one 
for the few-nucleon problem. These issues will be taken 
up in Sec. III. 

For more extended discussions of scales in the nuclear 
force problem and in nuclei, we refer to (10-13). 


C. Conventional approaches to the nuclear force problem 


Before discussing the application of the effective field 
theory approach to the nuclear force problem, let us make 
a few comments on the highly successful conventional 
approaches. First, we consider the two-nucleon case. 
Historically, meson field theory and dispersion relations 
have laid the foundations for the construction of a two- 
nucleon potential. All these approaches incorporate the 
long-range one-pion exchange as proposed by Yukawa in 
1935 (14) which nowadays is firmly rooted in QCD. Dis- 
persion relations can be used to construct the two-pion 
exchange contribution to the nuclear force as pioneered 
at Paris (15) and Stony Brook (16). For a review, see e.g. 
(17). In the 1990ties, the so-called high-precision poten- 
tials have been developed that fit the large basis of pp and 
np elastic scattering data with a y?/datum ~ 1. One of 
these is the so-called CD-Bonn potential (18) (which was 
developed at Moscow, Idaho). Besides one-pion, p and w 
vector-meson exchanges, it contains two scalar—isoscalar 
mesons in each partial wave up to angular momentum 


J = 5 with the mass and coupling constant of the second 
o fine-tuned in any partial wave. The hadronic vertices 
are regulated with form factors with cut-offs ranging from 
1.3 to 1.7 GeV. Similarly, in the Nijmegen I,II poten- 
tials one—pion exchange is supplemented by heavy boson 
exchanges with adjustable parameters which are fitted 
for all (low) partial waves separately (19). The Argonne 
V18 (AV18) potential starts from a very general oper- 
ator structure in coordinate space and has fit functions 
for all these various operators (20). While these various 
potentials give an accurate representation of the nucleon- 
nucleon phase shifts and of most deuteron properties, the 
situation becomes much less satisfactory when it comes 
to the much smaller but necessary three-nucleon forces. 
Such three-body forces are needed to describe the nuclear 
binding energies and levels, as most systematically shown 
by the Urbana-Argonne group (21). Systematic studies 
of the dynamics and reactions of systems with three or 
four-nucleons further sharpen the case for the necessity of 
including three-nucleon forces (3NFs), see e.g. (22). The 
archetype of a 3NF is due to Fujita and Miyazawa (FM) 
(23), who extended Yukawa’s meson exchange idea by 
sandwiching the pion-nucleon scattering amplitude be- 
tween nucleon lines, thus generating the 3NF of longest 
range. In fact, the work of Fujita and Miyazawa has been 
the seed for many meson-theoretical approaches to the 
three-nucleon force like the families of Tucson-Melbourne 
(24; 25), Brazilian (26) or Urbana-Illinois (27; 28) 3NFs. 
While the conventional approach as briefly outlined 
here as enjoyed many successes and is frequently used 
in e.g. nuclear structure and reaction calculations, it re- 
mains incomplete as there are certain deficiencies that 
can only be overcome based on EFT approaches. These 
are: (i) it is very difficult - if not impossible - to as- 
sign a trustworthy theoretical error, (ii) gauge and chi- 
ral symmetries are difficult to implement, (iii) none of 
the three-nucleon forces is consistent with the underly- 
ing nucleon-nucleon interaction models/approaches and 
(iv) the connection to QCD is not at all obvious. Still, 
as we will show later, there is a very natural connec- 
tion between these models and the forces derived from 
EFT by mapping the complicated physics of the short- 
distance part of any interaction at length scales ~ 1/M, 
to the tower of multi-fermion contact interactions that 
naturally arise in the EFT description (see Sec. II.B). 


D. Brief introduction to effective field theory 


Effective field theory (EFT) is a general approach to 
calculate the low-energy behavior of physical systems by 
exploiting a separation of scales in the system (for reviews 
see e.g. (29-31)). Its roots can be traced to the renormal- 
ization group (32) and the intuitive understanding of ul- 
traviolet divergences in quantum field theory (33). A suc- 
cinct formulation of the underlying principle was given 
by Weinberg (34): If one starts from the most general 
Lagrangian consistent with all symmetries of the under- 


FIG. 3 Expansion of a heavy-particle exchange diagram in 
terms of local light-particle operators. The solid and dashed 
lines denote light and heavy particles, respectively. The filled 
circle and square denote insertions with zero and two deriva- 
tives, in order. The ellipses stands for operators with more 
derivatives. 


lying interaction, one will get the most general S-matrix 
consistent with these symmetries. Together with a power 
counting scheme that specifies which terms are required 
at a desired accuracy leads to a predictive paradigm for 
a low-energy theory. The expansion is typically in pow- 
ers of a low-momentum scale Mjo. which can be the 
typical external momentum over a high-momentum scale 
Mhignh. However, what physical scales Mnignh and Miow 
are identified with depends on the considered system. In 
its most simple setting, consider a theory that is made 
of two particle species, the light and the heavy ones with 
Miow < Mnign. Consider now soft processes in which the 
energies and momenta are of the order of the light parti- 
cle mass (the so-called soft scale). Under such conditions, 
the short-distance physics related to the heavy particles 
can never be resolved. However, it can be represented 
by light-particle contact interactions with increasing di- 
mension (number of derivatives). Consider e.g. heavy 
particle exchange between light ones in the limit that 
Mhigh — © while keeping the ratio g/Mnign fixed, with g 
the light-heavy coupling constant. As depicted in Fig. 3, 
one can represent such exchange diagrams by a sum of 
local operators of the light fields with increasing number 
of derivatives. In a highly symbolic notation 


2 2 2 
g g gt 
——— = + Akers, (1.11) 
Mirign —t Mirign Mrrign 


with t the squared invariant momentum transfer. In 
many cases, the corresponding high-energy theory is not 
known. Still, the framework of EFT offers a predictive 
and systematic framework for performing calculations in 
the light particle sector. Denote by Q a typical energy 
or momentum of the order of Mjow and by A the hard 
scale where the EFT will break down. In many cases, 
this scale is set by the masses of the heavy particles not 
considered explicitely. In such a setting, any matrix ele- 
ment or Greens function admits an expansion in the small 
parameter Q/A (34) 


Q\" Q 
ma¥(k) 7s) 
where F is a function of order one (naturalness), ju a reg- 


ularization scale (related to the UV divergences appear- 
ing in the loop graphs) and the g; denotes a collection 


(1.12) 


of coupling constants, often called low-energy constants 
(LECs). These parameterize (encode) the unknown high- 
energy (short-distance) physics and must be determined 
by a fit to data (or can be directly calculated if the cor- 
responding high-energy theory is known/can be solved). 
The counting index v in general depends on the fields 
in the effective theory, the number of derivatives and 
the number of loops. This defines the so-called power 
counting which allows to categorize all contributions to 
any matrix element at a given order. It is important 
to stress that v must be bounded from below to define 
a sensible EFT. In QCD e.g. this is a consequence of 
the spontaneous breaking of its chiral symmetry. The 
contributions with the lowest possible value of v define 
the so-called leading order (LO) contribution, the first 
corrections with the second smallest allowed value of v 
the next-to-leading order (NLO) terms and so on. In 
contrast to more conventional perturbation theory, the 
small parameter is not a coupling constant (like, e.g., in 
Quantum Electrodynamcis) but rather one expands in 
small energies or momentum, where small refers to the 
hard scale A. The archetype of such a perturbative EFT 
is chiral perturbation theory that exploits the strictures 
of the spontaneous and explicit chiral symmetry break- 
ing in QCD (35; 36). Here, the light degrees of freedom 
are the pions, that are generated through the symmetry 
violation. Heavier particles like e.g. vector mesons only 
appear indirectly as they generate local four-pion interac- 
tions with four, six, ... derivatives. For a recent review, 
see Ref. (37). Of course, the pions also couple to heavy 
matter fields like e.g. nucleons, that can also be included 
in CHPT, as reviewed by Bernard (38). 

So far, we have made the implicit assumption of nat- 
uralness, which implies e.g. that the scattering length 
a is of natural size as e.g. in CHPT, where the scale 
is set by 1/Ay ~ 1GeV~' ~ 0.2fm. This also implies 
that there are no bound states close to the scattering 
thresholds. In many physical systems and of particular 
interest here, especially in the two-nucleon system, this is 
not the case, but one rather has to deal with unnaturally 
large scattering lengths (and also shallow bound states). 
To be specific, let us consider nucleon-nucleon scattering 
at very low energies in the 1So channel, cf. Eq. (1.10). 
For such low energies, even the pions can be considered 
heavy and are thus integrated out. To construct an EFT 
that is applicable for momenta p > 1/a, one must retain 
all terms ap ~ 1 in the scattering matrix. This requires 
a non-perturbative resummation and is most elegantly 
done in the power divergence scheme of Kaplan, Savage 
and Wise (39; 40). This amounts to summing the lead- 
ing four-nucleon contact term ~ Co(y'w)? to all orders 
in Co and matching the scale-dependent LEC Co to the 
scattering length. This leads to the T-matrix 


(1.13) 


where the expansion around the large scattering length 
is made explicit. All other effects, like e.g. effective 


range corrections, are treated perturbatively. This com- 
pact and elegant scheme is, however, not sufficient for 
discussing nuclear processes with momenta p > M,. We 
will come back to this topic when we give the explicit con- 
struction of the chiral nuclear EFT in sec. II.B. It is im- 
portant to stress that such EFTs with unnaturally large 
scattering length can exhibit universal phenomena that 
can be observed in physical systems which differ in their 
typical energy scale by many orders of magnitude, for a 
review see Braaten and Hammer (41). We also remark 
that there are many subtleties in constructing a proper 
EFT, but space forbids to discuss these here. Whenever 
appropriate and/or neccessary, we will mention these in 
the following sections and provide explicit references. 


E. First results from lattice QCD 


Lattice QCD (LQCD) is a promising tool to calculate 
hadron properties ab initio from the QCD Lagrangian on 
a discretized Euclidean space-time. This requires state- 
of-the-art high performance computers and refined algo- 
rithms to analyse the QCD partition function by Monte 
Carlo methods. Only recently soft- and hardware devel- 
opments have become available that allow for full QCD 
simulations at small enough quark masses (correspond- 
ing to pion masses below 300 MeV), large enough vol- 
umes (corresponding to spatial dimensions larger than 
2.5 fm) and sufficiently fine lattice spacing (a ~ 0.05 fm) 
so that the results are not heavily polluted by computa- 
tional artefacts and can really be connected to the phys- 
ical quark masses by sensible chiral extrapolations. 

For the nuclear force problem, there are two main de- 
velopments in LQCD to be reviewed here. These con- 
cern the extraction of hadron-hadron scattering lengths 
from unquenched simulations and the first attempts to 
construct a nuclear potential. These are groundbreaking 
studies, but clearly at present one has not yet achieved an 
accuracy to obtain high-precision predictions for nuclear 
properties. We look very much forward to the develop- 
ment of these approaches in the years to come. 

The first exploratory study of the nucleon-nucleon 
scattering lengths goes back to Fukugita et al. (42; 43) 
in the quenched approximation. They make use of an 
elegant formuala, frequently called the “Ltischer for- 
mula”, that relates the S-wave scattering length ag be- 
tween two hadrons h; and hg to the energy shift 6E = 
Ehihy — (m1 + mz) of the two-hadron state at zero rela- 
tive momentum confined in spatial box of size L°. It is 
given by (44-46) 


27 a0 ao aa =16 


with pw = myme2/(m, + m2) the reduced mass and c; = 
2.837297 and co = 6.375183. A generalization of this 
formalism was given by Beane et al.(47) utilizing methods 
developed for the so-called pionless nuclear EFT (EFT 


with contact interactions, for a review see e.g. (48)). It 
reads 


p cot do(p) = —S ((Lp/2n)?) , 
Ay 
S(n) = » prog oe (1.15) 


which gives the location of all energy eigenstates in the 
box. Here, do is the S-wave phase shift. The sum 
over all three-vectors of integers 7 is such that |7| < 
A; and the limit A; — oo is implicit. In the limit 
L > |ao|, Eq. (1.15) reduces to the Liischer formula, 
Eq. (1.14). On the other hand, for large scattering length, 
|(p cot do)~!| >> L, the energy of the lowest state is given 
by 
2 
Eo = a [dy + dz Lp cot 60 es | ; 
BL? 

with d, = 0.472895, dz = 0.0790234 and pcot do is evalu- 
ated at the energy E = 4n7d,/(uL). Within this frame- 
work, in Ref. (49) the first fully dynamical simulation 
of the neutron-proton scattering lengths was performed, 
with a lowest pion mass of 354 MeV. This mass is still 
too high to perform a precise chiral extrapolation to the 
physical pion mass, but this calculation clearly demon- 
strates the feasilbilty of this approach (see also sec. II.H). 
This scheme can also be extended to hyperon-nucleon in- 
teractions, see (50). A first signal for 7A and nX7 scat- 
tering was reported in Ref. (51). For a recent review on 
these activities of the NPLQCD collaboration, we refer 
the reader to (52). 

Another interesting development was initiated and car- 
ried out by Aoki, Hatsuda and Ishii (53). They have gen- 
eralized the two-center Bethe-Salpeter wavefunction ap- 
proach of the CP-PACS collaboration (54), which offers 
an alternative to Liischer’s formula, to the two-nucleon 
(NN) system. Given an interpolating field for the neu- 
tron and for the proton, the NN potential can be defined 
from the properly reduced 6-quark Bethe-Salpeter ampli- 
tude ¢(7). The resulting Lippmann-Schwinger equation 
defines a non-local potential for a given, fixed separation 
r = |r|. Performing a derivative expansion, the central 
potential V.(r) at a given energy E is extracted from 


(1.16) 


(1.17) 


Monte Carlo Simulations are then performed to gener- 
ate the 6-quark Bethe-Salpeter amplitude in a given spin 
and isospin state of the two-nucleon system on a large 
lattice, V = (4.4fm)* in the quenched approximation for 
pion masses between 380 and 730 MeV. Despite these ap- 
proximations, the resulting effective potential extracted 
using Eq. (1.17) shares the features of the phenomeno- 
logical NN potentials - a hard core (repulsion) at small 
separation surrounded by an attractive well at interme- 
diate and larger distances, see Fig. 4. Furthermore, the 
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FIG. 4 Effective potential in the 15> channel for three differ- 
ent quark masses in quenched LQCD after Ref. (55). The 
dashed line is the asymptotic OPEP for M, = 380MeV, 
m = 1.2GeV and g?y/(4r) = 14.0. We are grateful to Dr. 
N. Ishii for providing us with the data. 


asymptotic form of this potential has exactly the form of 
the OPE, provided one rescales the formula 


2m 


= = 2 .—M 
01:02 (=) e Mar 


with the pion and nucleon masses used in the simula- 
tions but keeping the pion-nucleon coupling at its phys- 
ical value, g2y/(4m) ~ 14.0. These interesting results 
have led to some enthusiastic appraisal, see e.g. (56). 
However, it is important to stress that the so-calculated 
potential is not unique, especially its properties at short 
distances, since it depends on the definition of the in- 
terpolating nucleon fields. Furthermore, the quenched 
approximation is known to have uncontrolled system- 
atic uncertainties as it does not even define a quantum 
field theory. In this context, the authors of Ref. (54) 
report on the numerical absence of the large-distance- 
dominating 7-exchange from the flavor-singlet hair-pin 
diagram. Still, one would like to see this promising cal- 
culation repeated with dynamical quarks of sufficiently 
small masses. In Ref. (57) this framework was used to 
study the =N interaction. Interestingly, the central po- 
tential of the p=° interaction looks very similar to the 
central np potential. It would be interesting to extend 
these calculations to other hyperon-nucleon channels and 
also study the effects of SU(3) symmetry breaking. We 
will come back to these issues in the context of an three- 
flavor chiral EFT in sec. II.F. For recent developments in 
this scheme concerning the calculation of the tensor force, 
the energy dependence of the NN potential and prelimi- 
nary results for full QCD (2+1 flavors), see the talks by 
Aoki, Ishii, and Nemura at the Lattice 2008 conference 
(58). 


F. Observables and not-so observable quantities 


There is an extensive literature, primarily from the 
sixties and seventies, on the role of off-shell physics in 
nuclear phenomena (see, e.g., Ref. (59) and references 
therein). This includes not only few-body systems (e.g., 
the triton) and nuclear matter, but interactions of two- 
body systems with external probes, such as nucleon- 
nucleon bremstrahlung and the electromagnetic form fac- 
tors of the deuteron. The implicit premise was that there 
is a true underlying potential governing the nucleon- 
nucleon force, so that its off-shell properties can be deter- 
mined. Indeed, the nuclear many-body problem has tra- 
ditionally been posed as finding approximate solutions to 
the many-particle Schrodinger equation, given a funda- 
mental two-body interaction that reproduces two-nucleon 
observables. 

In contrast, effective field theories are determined com- 
pletely by on-energy-shell information, up to a well- 
defined truncation error. In writing down the most gen- 
eral Lagrangian consistent with the symmetries of the un- 
derlying theory, many-body forces arise naturally. Even 
though they are usually suppressed at low energies, they 
enter at some order in the EFT expansion. These many- 
body forces have to be determined from many-body data. 
The key point is, however, that no off-energy-shell infor- 
mation is needed or experimentally accessible. A funda- 
mental theorem of quantum field theory states that phys- 
ical observables (or more precisely, S-matrix elements) 
are independent of the choice of interpolating fields that 
appear in a Lagrangian (60; 61). Equivalently, observ- 
ables are invariant under a change of field variables in a 
field theory Lagrangian (or Hamiltonian): 


Ve) > b(@) + nPldI, 


where P[w] is a local polynomial of the field 7 and 
its derivatives and 7 is an arbitray counting parame- 
ter. Newly generated contributions to observables have 
to cancel separately at each order in 7. This “equiva- 
lence theorem” holds for renormalized field theories. In 
an EFT, one exploits the invariance under field redefini- 
tions to eliminate redundant terms in the effective La- 
grangian and to choose the most convenient or efficient 
form for practical calculations (62-66). Since off-shell 
Green’s functions and the corresponding off-shell ampli- 
tudes do change under field redefinitions, one must con- 
clude that off-shell properties are unobservable. 

Several recent works have emphasized from a field 
theory point of view the impossibility of observing off- 
shell effects. In Refs. (67; 68), model calculations were 
used to illustrate how apparent determinations of the 
two-nucleon off-shell T-matrix in nucleon-nucleon brem- 
strahlung are illusory, since field redefinitions shift contri- 
butions between off-shell contributions and contact inter- 
actions. Similarly, it was shown in Ref. (69) that Comp- 
ton scattering on a pion cannot be used to extract infor- 
mation on the off-shell behavior of the pion form factor. 
The authors of Refs. (70; 71) emphasized the nonunique- 
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ness of chiral Lagrangians for three-nucleon forces and 
pion production. Field redefinitions lead to different off- 
shell forms that yield the same observables within a con- 
sistent power counting. In Ref. (72), an interaction pro- 
portional to the equation of motion is shown to have no 
observable consequence for the deuteron electromagnetic 
form factor, even though it contributes to the off-shell 
T-matrix. 

In systems with more than two nucleons, one can trade 
off-shell, two-body interactions for many-body forces. 
This explains how two-body interactions related by uni- 
tary transformations can predict different binding ener- 
gies for the triton (73) if many-body forces are not con- 
sistently included. These issues were discussed from the 
viewpoint of unitary transformations in Refs. (74) and 
(75). The extension to many-fermion systems in the ther- 
modynamic limit was considered in (76). The effects of 
field redefinitions were illustrated using the EFT for the 
dilute Fermi gas (77). If many-body interactions gener- 
ated by the field redefinitions are neglected, a Coester 
line similar to the one observed for nuclear matter (78) 
is generated. Moreover, the connection to more tradi- 
tional treatments using unitary transformations was elu- 
cidated. The question of whether occupation numbers 
and momentum distributions of nucleons in nuclei are 
observables was investigated in Ref. (79). Field redefini- 
tions lead to variations in the occupation numbers and 
momentum distributions that imply the answer is neg- 
ative. The natural size of the inherent ambiguity (or 
scheme dependence) in these quantities is determined by 
the applicability of the impulse approximation. Only if 
the impulse approximation is well justified, the ambiguity 
is small and these quantities are approximately scheme 
independent. This has important implications for the in- 
terpretation of (e,e’p) experiments with nuclei. Whether 
the stark difference in occupation numbers between non- 
relativistic and relativistic Brueckner calculations can be 
explained by this ambiguity is another interesting ques- 
tion (80). 


Il. EFT FOR FEW-NUCLEON SYSTEMS: 
FOUNDATIONS AND APPLICATIONS 


A. EFT with contact interactions and universal aspects 


In nuclear physics, there are a number of EFTs which 
are all useful for a certain range of systems (cf. Fig. 2). 
The simplest theories include only short range interac- 
tions and even integrate out the pions. At extremely low 
energies, Mpign is given by the NN scattering lengths 
and one can formulate a perturbative EFT in powers of of 
the typical momentum k divided by Mnign. Since the NN 
scattering lengths are large this theory has a very limited 
range of applicability. It is therefore useful to construct 
another EFT with short-range interactions that resums 
the interactions generating the large scattering length. 
This so-called pionless EFT can be understood as an ex- 


pansion around the limit of infinite scattering length or 
equivalently around threshold bound states. Its break- 
down scale is set by one-pion exchange, Mnhign ~ Mz, 
while Miow ~ 1/a ~ k. For momenta k of the order of 
the pion mass /,, pion exchange becomes a long-range 
interaction and has to be treated explicitly. This leads 
to the chiral EFT whose breakdown scale Mpjign is set 
by the chiral symmetry breaking scale A, and will be 
discussed in detail below. 

The pionless theory relies only on the large scattering 
length and is independent of the mechanism responsi- 
ble for it. It is very general and can be applied in sys- 
tems ranging from ultracold atoms to nuclear and parti- 
cle physics. It is therefore ideally suited to unravel uni- 
versal phenomena driven by the large scattering length 
such as limit cycle physics (81; 82) and the Efimov effect 
(83). For recent reviews of applications to the physics 
of ultracold atoms, see Refs. (41; 84). Here we consider 
applications of this theory in nuclear physics. 

The pionless EFT is designed to reproduce the well 
known effective range expansion. The leading order La- 
grangian can be written as: 
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where the dots represent higher-order terms suppressed 
by derivatives and more nucleon fields. The Pauli matri- 
ces 0; (Ta) operate in spin (isospin) space, respectively. 
The contact terms proportional to Ci (Cé) correspond to 
two-nucleon interactions in the 9S; (4S 9) NN channels. 
Their renormalized values are related to the correspond- 
ing large scattering lengths a; and a, in the spin-triplet 
and spin-singlet channels, respectively. The exact rela- 
tion, of course, depends on the renormalization scheme. 
Various schemes can be used, such as a momentum cutoff 
or dimensional regularization. Convenient schemes that 
have a manifest power counting at the level of individual 
diagrams are dimensional regularization with PDS sub- 
traction, where poles in 2 and 3 spatial dimensions are 
subtracted (40), or momentum subtractions schemes as 
in (85). However, a simple momentum cutoff can be used 
as well. 

Since the scattering lengths are set by the low- 
momentum scale a ~ 1/Miow, the leading contact in- 
teractions have to be resummed to all orders (40; 86). 
The nucleon-nucleon scattering amplitude in the 35, 
(*So) channels is obtained by summing the so-called bub- 
ble diagrams with the Cj (C§) interactions shown in 
Fig. 5. This summation gives the exact solution of the 
Lippmann-Schwinger equation for the C§ or C§ interac- 
tions. Higher order derivative terms which are not shown 
explicitly in Eq. (2.1) reproduce higher order terms in 
the effective range expansion. Since these terms are nat- 
ural and their size is set by Mnign, their contribution at 
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FIG. 5 The bubble diagrams with the contact interaction Cé 
or Co contributing to the two-nucleon scattering amplitude. 


low energies is suppressed by powers of Miow/Mnigh and 
can be treated in perturbation theory. The subleading 
correction is given by the effective range ro ~ 1/Mnigh 
and the corresponding diagrams are illustrated in Fig. 6. 
The renormalized S-wave scattering amplitude to next- 
to-leading order in a given channel then takes the form 


At 1 
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where k is the relative momentum of the nucleons and 
the dots indicate corrections of order (Miow/Mnign)? for 
typical momenta k ~ Miow. The pionless EFT becomes 
very useful in the two-nucleon sector when external cur- 
rents are considered and has been applied to a variety of 
electroweak processes. These calculations are reviewed in 
detail in Refs. (48; 87). More recently Christlmeier and 
GrieBhammer have calculated low-energy deuteron elec- 
trodisintegration in the framework of the pionless EFT 
(88). For the double differential cross sections of the 
d(e,e’) reaction at 0 = 180° excellent agreement was 
found with a recent experiment at S-DALINAC (89).#! 
The double-differential cross section for an incident elec- 
tron energy Ey = 27.8 MeV and 6 = 180° is shown in 
Fig. 7. The data were used to precisely map the M1 re- 
sponse which governs the reaction np — dy relevant to 
big-bang nucleosythesis. Finally, the reaction pp > ppr° 
near threshold was studied by Ando (90). 


We now proceed to the three-nucleon system. Here it 
is convenient (but not mandatory) to rewrite the theory 
using so-called “dimeron” auxilliary fields (91). We need 
two dimeron fields, one for each S-wave channel: (i) a 
field t; with spin (isospin) 1 (0) representing two nucleons 
interacting in the 3S; channel (the deuteron) and (ii) a 
field s, with spin (isospin) 0 (1) representing two nucleons 
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FIG. 6 Diagrams for the inclusion of higher order contact 
interactions. 


#1 However, there is a disagreement between theory and data for 
the small longitudinal-transverse interference contribution o,7 
that is currently not understood (88). 
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FIG. 7 Double-differential cross sections of the 7H(e, e’) re- 
action with errors (hatched bands) extracted from the ex- 
periment. The gray bands and dashed lines are calculations 
in pionless EFT and a potential model. Figure courtesy of 
H.W. GrieShammer. 


interacting in the 1.S9 channel (92): 
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where 7,7 are spin and a,b are isospin indices while g;, 
gs, At, As and G3 are the bare coupling constants. This 
Lagrangian goes beyond leading order and already in- 
cludes the effective range terms. The coupling constants 
gt, At, gs, As are matched to the scattering lengths aa, 
and effective ranges roq in the two channels (a = s,t). 
Alternatively, one can match to the position of the bound 
state/virtual state pole y,. in the T-matrix instead of the 
scattering length which often improves convergence (93). 
The two quantities are related through: 


1 
Vo =— (1 — 4/1 = 2700 /aa) ; 


T0a 


(2.4) 


where a = s,t. The term proportional to G3 consti- 
tutes a Wigner-SU(4) symmetric three-body interaction. 
It only contributes in the spin-doublet S-wave channel. 
When the auxilliary dimeron fields t; and s, are inte- 
grated out, an equivalent form containing only nucleon 
fields is obtained. At leading order when the effective 
range corrections are neglected, the spatial and time 
derivatives acting on the dimeron fields are omitted and 
the field is static. The coupling constants gq and Aa, 
a = s,t are then not independent and only the combina- 
tion g2/Aq enters in observables. This combination can 
then be matched to the scattering length or pole position. 

The simplest three-body process to consider is 
neutron-deuteron scattering below the breakup thresh- 


10 


° x 
+ al > OX 
FIG. 8 The integral equation for the boson-dimeron scatter- 


ing amplitude. The single (double) line indicates the boson 
(dimeron) propagator. 


old. In order to focus on the main aspects of renor- 
malization, we suppress all spin-isospin indices and com- 
plications from coupled channels in the three-nucleon 
problem. This corresponds to a system of three spin- 
less bosons with large scattering length. If the scattering 
length is positive, the bosons form a two-body bound 
state analog to the deuteron which we call dimeron. The 
leading order integral equation for boson-dimeron scat- 
tering is shown schematically in Fig. 8. For total orbital 
angular momentum L = 0, it takes the following form: 


16 
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Here, H determines the strength of the three-body force 
G3(A) = 2mH(A)/A? which enters already at leading or- 
der and A is a UV cutoff introduced to regularize the inte- 
gral equation. The magnitude of the incoming (outgoing) 
relative momenta is k (p) and E = 3k?/(4m) — 1/(ma?). 
The on-shell point corresponds to k = p and the phase 
shift can be obtained via k cot 6 = 1/T3(k,k; E)+ik. For 
H =0 and A > o, Eq. (2.5) reduces to the STM equa- 
tion first derived by Skorniakov and Ter-Martirosian (94). 
It is well known that the STM equation has no unique 
solution (95). The regularized equation has a unique so- 
lution for any given (finite) value of the ultraviolet cutoff 
A but the amplitude in the absence of the three-body 
force shows an oscillatory behavior on In A. Cutoff inde- 
pendence of the amplitude is restored by an appropriate 
“running” of H(A) which turns out to be a limit cycle 
(96; 97): 


cos[so In(A/A,.) + arctan so] 


AH A Ss 
( cos[so In(A/A,.) — arctan so] ’ 


(2.7) 
where A, is a dimensionful three-body parameter gener- 
ated by dimensional transmutation. Adjusting A, to a 
single three-body observable allows to determine all other 


low-energy properties of the three-body system. Note 
that the choice of the three-body parameter A, is not 
unique and there are other definitions more directly re- 
lated to experiment (41). Because H(A) in Eq. (2.7) 
vanishes for certain values of the cutoff A it is possible 
to eliminate the explicit three-body force from the equa- 
tions by working with a fixed cutoff that encodes the 
dependence on A,. This justifies tuning the cutoff A in 
the STM equation to reproduce a three-body datum and 
using the same cutoff to calculate other observables as 
suggested by Kharchenko (98). Equivalently, a subtrac- 
tion can be performed in the integral equation (99; 100). 
In any case, one three-body input parameter is needed for 
the calculation of observables. A comprehensive study of 
the range corrections to the three-boson spectrum was 
carried out in Ref. (101). The authors showed that all 
range corrections vanish in the unitary limit due to the 
discrete scale invariance. While the corrections propor- 
tional to r9/a vanish trivially, this includes also the cor- 
rections proportional to K,7rq where kK, = \/mB3 is the 
binding momentum of the Efimov state fixed by the cho- 
sen renormalization condition. Moreover, they have cal- 
culated the corrections to the Efimov spectrum for finite 
scattering length. The range corrections are negligible for 
the shallow states but become important for the deeper 
bound states. 


The integral equations for the three-nucleon problem 
derived from the Lagrangian (2.3) are a generalization 
of Eq. (2.5). (For their explicit form and derivation, see 
e.g. Ref. (102).) For S-wave nucleon-deuteron scatter- 
ing in the spin-quartet channel only the spin-1 dimeron 
field contributes. This integral equation has a unique so- 
lution for A — oo and there is no three-body force in 
the first few orders. The spin-quartet scattering phases 
can therefore be predicted to high precision from two- 
body data (103; 104). In the spin-doublet channel both 
dimeron fields as well as the the three-body force in the 
Lagrangian (2.3) contribute (92). This leads to a pair of 
coupled integral equations for the T-matrix. Thus, one 
needs a new parameter which is not determined in the 
2N system in order to fix the (leading) low-energy be- 
havior of the 3N system in this channel. The three-body 
parameter gives a natural explanation of universal corre- 
lations between different three-body observables such as 
the Phillips line, a correlation between the triton binding 
energy and the spin-doublet neutron-deuteron scattering 
length (105). These correlations are purely driven by the 
large scattering length independent of the mechanism re- 
sponsible for it. As a consequence, they occur in atomic 
systems such as *He atoms as well (41). 


Higher-order corrections to the amplitude including 
the ones due to 2N effective range terms can be included 
perturbatively. This was first done at NLO for the scat- 
tering length and triton binding energy in (106) and for 
the energy dependence of the phase shifts in (99). In 
Refs. (102; 107), it was demonstrated that it is conve- 
nient to iterate certain higher order range terms in or- 
der to extend the calculation to N2LO. Here, also a sub- 
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FIG. 9 Phase shifts for neutron-deuteron scattering below 
the deuteron breakup at LO (dash-dotted line), NLO (dashed 
line), and N?LO (solid line). The filled squares and circles 
give the results of a phase shift analysis and a calculation 
using AV18 and the Urbana IX three-body force, respectively. 
Figure courtesy of L. Platter. 


leading three-body force was included as required by di- 
mensional analysis. More recently, Platter and Phillips 
showed using the subtractive renormalization that the 
leading three-body force is sufficient to achieve cutoff in- 
dependence up to N*LO in the expansion in Mjow/Mnign 
(108). The results for the spin-doublet neutron-deuteron 
scattering phase shift at LO (92), NLO (99), and N?LO 
(109) are shown in Fig. 9. There is excellent agreement 
with the available phase shift analysis and a calculation 
using a phenomenological NN interaction. Whether there 
is a suppression of the subleading three-body force or 
simply a correlation between the leading and subleading 
contributions is not fully understood. The extension to 
3N channels with higher orbital angular momentum is 
straightforward (110) and three-body forces do not ap- 
pear until very high orders. A general counting scheme 
for three-body forces based on the asymptotic behavior of 
the solutions of the leading order STM equation was pro- 
posed in (111). A complementary approach to the few- 
nucleon problem is given by the renormalization group 
where the power counting is determined from the scaling 
of operators under the renormalization group transforma- 
tion (32). This method leads to consistent results for the 
power counting (112-114). Universal low-energy prop- 
erties of few-body systems with short-range interactions 
and large two-body scattering length were reviewed in 
(41). (See also (115) for an early work on this subject.) 
Three-body calculations with external currents are still 
in their infancy. However, a few exploratory calculations 
have been carried out. Universal properties of the tri- 
ton charge form factor were investigated in Ref. (116) 
and neutron-deuteron radiative capture was calculated 
in Refs. (117; 118). This opens the possibility to carry 
out accurate calculations of electroweak reactions at very 
low energies for astrophysical processes. 


The pionless approach has also been extended to the 
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FIG. 10 The Tjon line correlation as predicted by the pionless 
theory. The grey circles and triangles show various calcula- 
tions using phenomenological potentials (121). The squares 
show the results of chiral EFT at NLO for different cutoffs 
while the diamond gives the N*LO result (122; 123). The 
cross shows the experimental point. 


four-body sector (119; 120). In order to be able to ap- 
ply the Yakubovsky equations, an equivalent effective 
quantum mechanics formulation was used. The study 
of the cutoff dependence of the four-body binding en- 
ergies revealed that no four-body parameter is required 
for renormalization at leading order. As a consequence, 
there are universal correlations in the four-body sector 
which are also driven by the large scattering length. The 
best known example is the Tjon line: a correlation be- 
tween the triton and alpha particle binding energies, By 
and By, respectively. Of course, higher order correc- 
tions break the exact correlation and generate a band. 
In Fig. 10, we show this band together with some some 
calculations using phenomenological potentials (121) and 
achiral EFT potential with explicit pions (122; 123). All 
calculations with interactions that give a large scatter- 
ing length must lie within the band. Different short- 
distance physics and/or cutoff dependence should only 
move the results along the band. This can for example 
be observed in the NLO results with the chiral potential 
indicated by the squares in Fig. 10 or in the few-body 
calculations with the low-momentum NN potential Viow x 
carried out in Ref. (124). The Viow x potential is obtained 
from phenomenological NN interactions by intergrating 
out high-momentum modes above a cutoff A but leav- 
ing two-body observables (such as the large scattering 
lengths) unchanged. The results of few-body calculations 
with Viowk are not independent of A but lie all close to 
the Tjon line (cf. Fig. 2 in Ref. (124)). The studies of 
the four-body system in the pionless theory were exten- 
dend further in Ref. (125). Here the dependence of the 
four-body bound state spectrum on the two-body scat- 
tering length was investigated in detail and summarized 
in a generalized Efimov plot for the four-body spectrum. 

The question of whether a four-body parameter has 
to enter at leading order was reanalyzed by Yamashita 
et al. (126). Within the renormalized zero-range model, 
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they found a strong sensitivity of the deepest four-body 
energy to a four-body subtraction constant in their equa- 
tions. They motivated this observation from a general 
model-space reduction of a realistic two-body interaction 
close to a Feshbach resonance. The results of Ref. (119) 
for the 4He tetramer that include a four-body parame- 
ter were also reproduced. Yamashita et al. concluded 
that a four-body parameter should generally enter at 
leading order. They argued that four-body systems of 
4He atoms and nucleons (where this sensitivity is ab- 
sent (119; 120; 124)) are special because repulsive in- 
teractions strongly reduce the probability to have four 
particles close together. However, the renormalization of 
the four-body problem was not explicitly verified in their 
calculation. Another drawback of their analysis is the fo- 
cus on the deepest four-body state only. Therefore, their 
findings could be an artefact of their particular regular- 
ization scheme. Another recent study by von Stecher and 
collaborators (127) confirmed the absence of a four-body 
parameter for shallow states while some sensitivity was 
found for the deepest four-body state. 

The pionless theory has also been extended to more 
than four particles using it within the no-core shell model 
approach. Here the expansion in a truncated harmonic 
oscillator basis is used as the ultraviolet regulator of the 
EFT. The effective interaction is determined directly in 
the model space, where an exact diagonalization in a 
complete many-body basis is performed. In Ref. (128), 
the 0+ excited state of *He and the °Li ground state 
were calculated using the deuteron, triton, alpha particle 
ground states as input. The first (0T;0) excited state in 
4He is calculated within 10% of the experimental value, 
while the °Li ground state comes out at about 70% of 
the experimental value in agreement with the 30 % error 
expected for the leading order approximation. These re- 
sults are promising and should be improved if range cor- 
rections are included. Finally, the spectrum of trapped 
three- and four-fermion systems was calculated using the 
same method (129). In this case the harmonic potential 
is physical and not simply used as an ultraviolet regula- 
tor. 


B. Chiral EFT for few nucleons: foundations 


The extension of the previously discussed EFT with 
contact interactions to higher energies requires the in- 
clusion of pions as explicit degrees of freedom. The in- 
teraction between pions and nucleons can be described 
in a systematic way using chiral perturbation theory. In 
contrast, the interaction between the nucleons is strong 
and leads to nonperturbative phenomena at low energy 
such as e.g. shallow-lying bound states. This break- 
down of perturbation theory can be linked to the fact 
that the interaction between the nucleons is not sup- 
pressed in the chiral limit contrary to the pion and pion- 
nucleon interactions. Moreover, an additional enhance- 
ment occurs for Feynman diagrams involving two and 


Ase 


FIG. 11 Representation of the two-pion exchange Feynman 
diagram in terms of time-ordered graphs. Solid and dashed 
lines represent nucleons and pions, respectively. 


more nucleons due to the appearance of the so-called 
pinch singularities in the limit of the infinite nucleon 
mass. Although such infrared singularities disappear if 
one keeps the nucleon mass at its physical value, they 
do generate large enhancement factors which destroy the 
chiral power counting. This can be more easily under- 
stood utilizing the language of time-ordered perturbation 
theory. Consider, for example, the two-pion exchange 
box diagram shown in Fig. 11. While all intermediate 
states in the first two time-ordered graphs, often referred 
to as irreducible, involve at least one virtual pion and 
thus lead to energy denimonators of the expected size, 
E— E, ~ M,, the remaining reducible diagrams in- 
volve an intermediate state with nucleons only which pro- 
duces unnaturally small energy denominators of the order 
E-E; ~ M2/m« M,. Clearly, the enhanced reducible 
time-ordered diagrams are nothing but the iterations of 
the Lippmann-Schwinger equation with the kernel which 
contains all possible irreducible diagrams and defines the 
nuclear Hamiltonian. It is free from infrared enhance- 
ment factors and can be worked out systematically using 
the machinery of chiral perturbation theory as suggested 
in Weinberg’s seminal work (130; 131).%?. This natural 
reduction to the quantum mechanical A-body problem 
is a welcome feature for practical calculations as it al- 
lows to apply various existing few-body techniques such 
as e.g. the Faddeev- Yakubovsky scheme, the no-core shell 
model, Green’s function Monte Carlo and hyperspherical 
harmonics methods. On the other hand, the framework 
offers a systematic and perturbative scheme to derive nu- 
clear forces and current operators in harmony with the 
chiral symmetry of QCD. The expansion parameter is 
given by the ratio Q/A where Q is the soft scale asso- 
ciated with the pion mass and/or external nucleon mo- 
menta and A is the pertinent hard scale. For a given 
connected irreducible diagram with N nucleons, L pion 
loops and V; vertices of type i, the power v of the soft 
scale Q which determines its importance can be obtained 
based on naive dimensional analysis (i.e. assuming clas- 


#2 An alternative framework based on the perturbative treatment of 
the pion exchange contributions has been introduced by Kaplan, 
Savage and Wise (39; 40), see (48; 87) for review articles. As 
shown in Refs. (132-134), the perturbative inclusion of pions 
does not allow to significantly increase the applicability range 
of the theory as compared to pionless EFT. For yet different 
proposals to include pions in EFT for the nucleon-nucleon system 
see Refs. (135-137). 
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sical scaling dimensions for various operators in the ef- 
fective Lagrangian): 


1 
y=-44+2N+2L4+) ViAi, A= dit5ni—2. (2.8) 


Here, n; is the number of nucleon field operators and 
d; the number of derivatives and/or insertions of M,. 
The spontaneously broken chiral symmetry of QCD guar- 
antees A; > 0. As a consequence, the chiral dimen- 
sion v is bounded from below, and only a finite num- 
ber of diagrams contribute at a given order. In addition, 
Eq. (2.8) provides a natural explanation to the domi- 
nance of the two-nucleon interactions and the hierarchy 
of nuclear forces observed in nuclear physics. In partic- 
ular, it implies that two-, three-, and four-nucleon forces 
start to contribute at orders v = 0, 2 and 4, respectively. 
Notice that as argued in Ref. (131), the nucleon mass m 
should be counted as Q/m ~ Q?/A? (which implies that 
m > A) in order to maintain consistency with the ap- 
pearance of shallow-lying bound states.7*? Notice further 
that according to this counting rule, the momentum scale 
associated with the real pion production is treated as the 
hard scale, mM, ~ A, and needs not be explicitly kept 
track of, see also (138) and section II.E for a related dis- 
cussion. Clearly, such a framework is only applicable at 
energies well below the pion production threshold. We 
also emphasize that the validity of the naive dimensional 
scaling rules for few-nucleon contact operators has been 
questioned in Refs. (139; 140). We will come back to this 
issue in section II.C. 

Before discussing the chiral expansion of the nuclear 
forces it is important to clarify the relation between the 
underlying chiral Lagrangian for pions and nucleons and 
the nuclear Hamiltonian we are finally interested in. The 
derivation of the nuclear potentials from field theory is an 
old and extensively studied problem in nuclear physics. 
Different approaches have been developed in the fifties 
of the last century in the context of the so-called me- 
son theory of nuclear forces, see e.g. the review article 
(141). In the modern framework of chiral EFT, the most 
frequently used methods besides the already mentioned 
time-ordered perturbation theory are the ones based on 
S-matrix and the unitary transformation. In the former 
scheme, the nuclear potential is defined through match- 
ing the amplitude to the iterated Lippmann-Schwinger 
equation (142). In the second approach, the potential is 
obtained by applying an appropriately choosen unitary 
transformation to the underlying pion-nucleon Hamilto- 
nian which eliminates the coupling between the purely 
nucleonic Fock space states and the ones which contain 
pions, see (143) for more details. We stress that both 
methods lead to energy-independent interactions as op- 
posed by the ones obtained in time-ordered perturbation 


#3 This statement only applies for the power counting based on 
naive dimensional analysis. 


theory. The energy independence of the potential is a 
welcome feature which enables applications to three- and 
more-nucleon systems. 

We are now in the position to discuss the structure of 
the nuclear force at lowest orders of the chiral expansion. 
The leading-order (LO) contribution results, according 
to Eq. (2.8), from two-nucleon tree diagrams constructed 
from the Lagrangian of lowest dimension A; = 0, £L, 
which has the following form in the heavy-baryon formu- 
lation (144; 145): 


F? 
m & 
£Q = N(iv-D+ Ga u-S)N, 


= 
I 


(VAUV UT + x4), (2.9) 


B 
2S 
Se 

lI 


~5Cs(NN)(WN) + 2Cr(NSN)- (NSN), 


where N, v, and S,, = (1/2)iyso,,v" denote the large 
component of the nucleon field, the nucleon four-velocity 
and the covariant spin vector, respectively. The brackets 


(...) denote traces in the flavor space while F' and GA 
refer to the chiral-limit values of the pion decay and the 
nucleon axial vector coupling constants. The low-energy 
constants (LECs) C's and Cr determine the strength of 
the leading NN short-range interaction. Further, the uni- 
tary 2x2 matrix U(m) = u?(z) in the flavor space collects 
the pion fields, 


Cs oe ee eee ee 


* ae (2.10) 


where 7; denotes the isospin Pauli matrix. The covariant 
derivatives of the nucleon and pion fields are defined via 
Dy = On + (ul, Opul/2 and u, = i(utd,u — ud,ut). The 
quantity y+ = ulyul +uxtu with y = 2BM involves the 
explicit chiral symmetry breaking due to the finite light 
quark masses, M = diag(mu, ma). The constant B is 
related to the value of the scalar quark condensate in the 
chiral limit, (0|au|0) = —F?B, and relates the pion mass 
M,, to the quark mass mg via M? = 2Bm,+O(m). For 
more details on the notation and the complete expres- 
sions for the pion-nucleon Lagrangian including up to 
four derivatives/M,-insertions the reader is referred to 
(146). Expanding the effective Lagrangian in Eqs. (2.9) 
in powers of the pion fields one can easily verify that 
the only possible connected two-nucleon tree diagrams 
are the one-pion exchange and the contact one, see the 
first line in Fig. 12, yielding the following potential in the 
two-nucleon center-of-mass system (CMS): 
2 61° GG2°¢ 
VO) = Srieerear ca pp Ce kOpe so), OAL) 
where the superscript of V2n denotes the chiral order v, 
o; are the Pauli spin matrices, g = p’ — p is the nucleon 
momentum transfer and p (p'’) refers to initial (final) 
nucleon momenta in the CMS. Further, F, = 92.4 MeV 
and g4 = 1.267 denote the pion decay and the nucleon 


axial coupling constants, respectively. 
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FIG. 12 Chiral expansion of the two-nucleon force up to 
N°LO. Solid dots, filled circles, squares and diamonds denote 
vertices with A; = 0, 1, 2 and 4, respectively. Only irre- 
ducible contributions of the diagrams are taken in to account 
as explained in the text. 


The first corrections to the LO result are suppressed 
by two powers of the low-momentum scale. The ab- 
sence of the contributions at order vy = 1 can be traced 
back to parity conservation which forbids (NN)(NN) 
vertices with one spatial derivative and 7NWN vertices 
with two derivatives (i. e. A; = 1). The next-to-leading- 
order (NLO) contributions to the 2NF therefore result 
from tree diagrams with one insertion of the A; = 2- 
interaction and one-loop diagrams constructed from the 
lowest-order vertices, see Fig. 12. The relevant terms in 
the effective Lagrangian read (147) 


l l 
Ona Ba ee nqyt 
LE) = Soa)? + (2 UV!) Ox4) 
+ 2x'Uxty + xUtyUt) — 4oxty) - (x-)?) 


eae) 


2 1 1 
LS = x(Sw -D)? ~~ D-D+digS- u(x+) 


2m 2m 


+ idigS*"[Dh, x] + . ai 


co), = —C;((NDN)-(NDN) + ((DN)N)- ((DN)N)) 
— 2(C, + C2)(NDN)-((DN)N) 


— Co(NN)-((D?N)N+ND?N)+..., 


(2.12) 


where l;, d; and C; denote further LECs and m is the 
nucleon mass in the chiral limit. The ellipses in the pion 
and pion-nucleon Lagrangians refer to terms which do 
not contribute to the nuclear force at. NLO. In the case 
of the nucleon-nucleon Lagrangian co. only a few terms 
are given explicitly. The complete reparametrization- 
invariant set of terms can be found in (148). The NLO 
contributions to the two-nucleon potential have been first 
considered in (149; 150) utilizing the framework of time- 
ordered perturbation theory. The corresponding energy- 
independent expressions have been worked out in (151) 
using the method described in (152) and then re-derived 
in (142) using an S-matrix-based approach and, inde- 
pendently, in (143; 153) based on the method of uni- 
tary transformation. The one-pion (17) exchange dia- 
grams at NLO do not produce any new momentum de- 
pendence. Apart from renormalization of various LECs 
in Eq. (2.11), one obtains the leading contribution to the 
Goldberger-Treiman discrepancy (154) 


Grn _ ga 2M; 


dig tas. 
ne Ee 


(2.13) 


where the ellipses refer to higher-order terms. Similarly, 
loop diagrams involving NN short-range interactions only 
lead to (M,-dependent) shifts in the LO contact terms. 
The remaining contributions to the 2NF due to higher- 
order contact interactions and two-pion exchange have 
the form: 


T1:T2 A 
Vi) = ~ ETE 1M a) (4Mz608 — 494-1) 
4894,M3 
2 4 2 Aen 
9364 — We, 1) 
394 he ae 2 
— St EN) (A1- P62 G- G1- H2g") 


+ C1 7? + Cok? + (C3 G2 + Cab’) B41 - Se 
ae pe A aS eo ee 
+ is 5 (F1 +02) FX K+ CoG- 1G: Fe 


> > 


+ Crk-ak- G2, 


(2.14) 


where g = |¢| and the LECs C; can be written as lin- 
ear combinations of C; in Eq. (2.12). The loop function 


LA(q) is defined in the spectral function regularization 
(SFR) (155; 156) as 


= . Aw? 252 1 9A 
Aw) = 6th —om,) 2. metre + 2haws 
2q 4M2(A? + q?) 
(2.15) 
where we have introduced the following abbreviations: 


w = \/4M2 + @? and s = ,/A? — 4M2?. Here, A denotes 
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FIG. 13 Chiral expansion of the three-nucleon force up to 
N°LO. Diagrams in the first line (NLO) yield vanishing con- 
tributions to the 3NF if one uses energy-independent for- 
mulations as explained in the text. The five topologies 
at N?LO involve the two-pion exchange, one-pion-two-pion- 
exchange, ring, contact-one-pion exchange and contact-two- 
pion-exchange diagrams in order. Shaded blobs represent 
the corresponding amplitudes. For remaining notation see 
Fig. 12. 


the ultraviolet cutoff in the mass spectrum of the two- 
pion-exchange potential. If dimensional regularization 
(DR) is employed, the expression for the loop function 
simplifies to 


(2.16) 


In addition to the two-nucleon contributions, at NLO 
one also needs to consider three-nucleon diagrams shown 
in the first line of Fig. 13. The first diagram does not in- 
volve reducible topologies and, therefore, can be dealt 
with using the Feynman graph technique. It is then 
easy to verify that its contribution is shifted to higher 
orders due to the additional suppression by the factor 
of 1/m caused by the appearance of time derivative at 
the leading-order mmNN vertex, the so-called Weinberg- 
Tomozawa vertex. The two remeining diagrams have 


been considered by Weinberg (130; 131) and later by Or- 
donez and van Kolck (157) using the energy-dependent 
formulation based on time-ordered perturbation theory. 
In this approach, it was shown that the resulting 3NF 
cancels exactly (at the order one is working) against the 
recoil correction to the 2NF when the latter is iterated 
in the dynamical equation. In energy-independent ap- 
proaches such as e.g. the method of unitary transforma- 
tion which are employed in most of the existing few- 
nucleon calculations one observes that the irreducible 
contributions from the last two diagrams in the first line 
of Fig. 13 are suppressed by the factor 1/m and thus 
occur at higher orders (148), see also (158; 159). Conse- 
quently, there is no 3NF at NLO in the chiral expansion. 

The contributions at next-to-next-to-leading order 
(N?LO) involve one-loop diagrams with one insertion of 
the subleading vertices of dimension A; = 1, see Fig. 12. 
The corresponding Lagrangians read: 


LR = N(cr (x4) tea (vu)? +esu-u 
+ ca [S", S”]uptty + es (X4))N , 
Do- = 
Lenn = 5(NN)(NS-uN), 

where ¥+ = x+ — (yx+)/2 and D, c; are the LECs. The 
1m-exchange loop diagram again only lead to renormal- 
ization of the corresponding LECs. Similarly, the con- 
tribution from the last diagram which involves the two- 
nucleon contact interaction can be absorbed into a re- 
definition of the LECs C's,r and C; in Eqs. (2.11), (2.14) 
(provided one is not interested in the quark mass depen- 
dence of the nuclear force). Further, the football diagram 
yields vanishing contribution due to the antisymmetric 
(with respect to pion isospin quantum numbers) nature 
of the Weinberg-Tomozawa vertex. Thus, the only non- 
vanishing contribution at this order results from the 27- 
exchange triangle diagram: 


(2.17) 


3 39% - 2 
vO = = len (2M?2(2c1 — c3) — cs) (2M? +g”) 
x AA(g) — BAO Ss (4M2 + q?)A*(q) 
q 320 FA 1 Te eS | q 
x (41 -@G2-¢-—G? GG), (2.18) 
where the loop function AA(q) is given by 
: , 1 A-2M, 
A*(q) = 6(A — 2M,,) — arctan dai) . (2.19) 
2q gq? + 2AM, 


In DR, the expression for A(q) takes the following simple 
form: 


A(q) = lim AX(q) vee arctan 


; 5 5 (2.20) 


Notice that the triangle diagram also generates short- 
range contributions which may be absorbed into re- 
definition of contact interactions. The isoscalar central 
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contribution proportional to the LEC cg is attractive 
and very strong. It is by far the strongest two-pion 
exchange contribution and reaches a few 10s of MeV 
(depending on the choice of regularization) at internu- 
cleon distances of the order r ~ M1. The origin of 
the unnaturally strong subleading 27-exchange contribu- 
tions can be traced back to the (numerically) large val- 
ues of the LECs c3,4 and is well understood in terms of 
resonance exchange related to A excitation (160). We 
will come back to this issue in section II.D where the 
chiral EFT formulation with explicit A degrees of free- 
dom will be discussed. The central 27-exchange potential 
was also calculated by Robilotta (161) using the infrared- 
reglarized version of chiral EFT which enables to sum up 
a certain class of relativistic corrections (162). He found 
that the results in the heavy-baryon limit overestimate 
the ones obtained using infrared regularization by about 
25%, see also (163) for a related discussion. Last but 
not least, the chiral 27-exchange potential up to N?LO 
has been tested in the Nijmegen partial wave analysis 
(PWA) of both proton-proton and neutron-proton data 
(164; 165) where also an attempt has been done to de- 
termine the values of the LECs c3 4. As demonstrated 
in these studies, the representation of the (strong) long- 
range interaction based on the combination of the I7- 
and the chiral 27-exchange potentials rather than on the 
pure 1z-exchnge potential allows to considerably reduce 
the number of phenomenological parameters entering the 
energy-dependent boundary conditions which are needed 
to parametrize the missing short- and medium-range in- 
teractions. Also the extracted values of the LECs c3,4 
agree reasonably well with various determinations in the 
pion-nucleon system. These studies provide a beautiful 
confirmation of the important role of the 27-exchange 
potential in nucleon-nucleon scattering observables, see, 
however, Ref. (166) for a criticism. For a similar work uti- 
lizing the distorted-wave methods the reader is referred 
to (167; 168). 

The first nonvanishing contributions to the 3NF also 
show up at N?LO and arise from tree diagrams shown in 
Fig. 13 which involve a single insertion of the subleading 
vertices £™) in Eq. (2.17) and 

1. = - 
LO = —5(NN)(NTN)- (NrN). (2.21) 
where F is a low-energy constant. 
3NF expression read: 


The corresponding 
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where the subscripts refer to the nucleon labels and 
dG; = pi' — pi, with p;’ and p; being the final and initial 


momenta of the nucleon 7. The expressions in Eq. (2.22) 
correspond to a particular choice of nucleon labels. The 
full expression for the 3NF results by taking into account 
all possible permutations of the nucleons (for three nu- 
cleons there are altogether six permutations), i.e.: 


ya = V3n + all permutations. (2.23) 


We further emphasize that the expressions for the 3NF 
given in Refs. (157; 169) contain one redundant 17- 
exchange and two redundant contact interactions. As 
shown in (123), only one independent linear combina- 
tion contributes in each case if one considers matrix el- 
ements between antisymmetrized few-nucleon states, see 
also (92) for a related discussion. 

We now turn to next-to-next-to-next-to-leading order 
(NLO) and discuss first the corrections to the 2NF. As 
follows from Eq. (2.8), one has to account for contribu- 
tions from tree diagrams with one insertion from £™ or 
two insertions from £), one-loop diagrams with one in- 
sertion from £°) or two insertions from £“ as well as 
two-loop graphs constructed from the lowest-order ver- 
tices, see Fig. 12. Apart from renormalization of various 
LECs, the 1m-exchange potential receives at this order (in 
the scheme based on the counting m ~ O(A?/M,,)) the 
first relativistic corrections proportional to m~?. These 
are scheme-dependent and have to be chosen consistently 
with the 1/m-corrections to the 27-exchange potential 
and the relativistic extension of the dynamical equation, 
see (170) for a comprehensive discussion. The two-pion 
exchange contributions at N?LO were worked out by 
Kaiser (171) based on the one-loop representation of the 
aN scattering amplitude. We refrain from giving here 
the rather involved expressions for the sub-subleading 27- 
exchange potential and refer to the original work (171) 
where the results are given in terms of the corresponding 
spectral functions. For certain classes of contributions, 
the integrals over the two-pion exchange spectrum could 
be performed analytically and are given in Ref. (172). 
Notice further that the subleading (i.e. the ones propor- 
tional to m~?) relativistic corrections of the 27-exchange 
range have also been worked out by Kaiser (173). In the 
counting scheme with m ~ O(A?/M,,), these terms, how- 
ever, would only appear at N°LO. It should also be em- 
phasized that the N°LO contributions to the 27-exchange 
potential were worked out in the covariant version of chi- 
ral EFT (more precisely, using the formulation by Becher 
and Leutwyler (162)) by Higa et al. (174-176). 

37-exchange contributions also appear at this order 
in the chiral expansion and have been worked out in 
Refs. (177; 178), see also Ref. (179) for a related work. 
The resulting potentials turn out to be rather weak. For 
example, the strongest contribution is of the isoscalar 
spin-spin type (i.e. proportional to @; - 2) and about 10 
times weaker than the corresponding 27-exchange contri- 
bution at the same order at relative distances r ~ M71. 
It should, however, be emphasized that the subleading 
37-exchange contributions at N4*LO are larger in size 
(180) which, again, can be traced back to the large val- 
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ues of the LECs c;. Finally, the last type of the 2NF 
corrections at this order results from diagrams involving 
contact interactions. The most general polynomial (in 
momenta) representation of the short-range part of the 
potential involves, apart from the two leading and seven 
subleading terms given in Eqs. (2.11) and (2.14) fifteen 
new contact interactions (in the isospin invariant sector) 
yielding in total 24 LECs to be determined from nucleon- 
nucleon data. 

The 3NF contributions at N°LO feed into five different 
topologies, see Fig. 13, and are currently being worked 
out. Currently, the expressions for the first three topolo- 
gies which do not involve short-range contact interac- 
tions are available. The one-loop corrections to the 27- 
exchange diagrams can, to a large extent, be accounted 
for by a finite shift c; — ¢ = c; + dc; of the LECs c; 
(181; 182) 


94 M, = gs M;, 
64n Fe?’ 


6c1 = (2.24) 
Numerically, these corrections are of the order of 20% of 
the corresponding LECs and are consistent with the dif- 
ference in values of c; between the order-Q? and Q? deter- 
minations from the pion-nucleon system, see (160; 183- 
185). The only 27-exchange contribution that cannot 
be cast into redefinition of the LECs c; arises from the 
diagram which involves pions interacting in flight, see 
(181; 182) for the explicit expression. We also empha- 
size that there are no 27-exchange contributions from 
tree diagrams with one insertion from £) in Eq. (2.12) 
(except for the relativistic corrections). This is because 
diagrams involving subleading 7NN interaction do not 
yield any irreducible contributions while the ones with 
the aaNN-vertices of dimension v = 2 involve at least 
one time derivative and are, therefore, suppressed by a 
factor of 1/m. This observation is consistent with the ab- 
sence of logarithmic ultraviolet divergences in the loop di- 
agrams. In this context, it should be emphasized that the 
requirement of renormalizability at N°LO (and, presum- 
ably, also at higher orders) was found to impose strong 
constraints on the unitary ambiguity in the form of the 
resulting nuclear potentials. This issue is discussed ex- 
tensively in (186) and may remind one of the recent find- 
ings in the context of large-N. QCD (187-189) where it 
was shown that the multiple-meson-exchange potential 
derived in the energy-dependent formulation is incon- 
sistent with large-Nc counting rules. The consistency 
could be maintained using a different (but equivalent 
up to the considered order) form of the potential based 
on the energy-independent formalism, see Ref. (189) for 
more details. The contributions from the two-pion-one- 
pion exchange and ring diagrams are given explicitly in 
Ref. (181) where expressions are shown both in momen- 
tum and coordinate spaces. Especially in the case of ring 
diagrams where loop integrals involve two independent 
external momenta and, therefore, yield rather involved 
expressions in momentum space. It is advantageous to 
switch to coordinate space where a much more compact 
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FIG. 14 Diagrams contributing to the four-nucleon force at 
N®LO. For notation see Fig. 12. 


representation emerges. Notice further that ring dia- 
grams were already studied in the pioneering work (190). 
The calculation of the last two topologies involving the 
leading contact interactions is in progress. Last but not 
least, one should also take into account the leading rel- 
ativistic 1/m-corrections to the NLO three-nucleon dia- 
grams, see the first line in Fig. 13. Again, these contri- 
butions are scheme-dependent and should be chosen con- 
sistently with the relativistic corrections to the 2NF and 
the form of the dynamical equation. The 1/m-corrections 
to the 27-exchange 3NF have already been worked out 
long time ago by Coon and Friar and are given in the 
most general form in Ref. (158). Notice further that at 
this order, one needs to account for the dependence of 
the 2NF on the total momentum of the 2N system (ef- 
fects due to drift of the CMS of a two-body subsystem). 
Such boosted 2N operators may, in fact, also be viewed 
as 3N operators. In the context of chiral EFT, this kind 
of corrections is discussed in (191). 


The last type of N?LO contributions arises from four- 
nucleon tree diagrams constructed from the lowest-order 
vertices, see in Fig. 14, which have been evaluated re- 
cently using the method of unitary transformation (186; 
192). Notice that the first two diagrams in the second 
line were already discussed long ago, see e.g. (193; 194). 
Furthermore, it should be emphasized that disconnected 
diagrams calculated e.g. in (169) using time-ordered per- 
turbation theory do not contribute to the nuclear force in 
the method of unitary transformation. It has been con- 
jectured in Ref. (194) that 4N diagrams which involve 
reducible topologies do not generate irreducible pieces in 
the amplitude and thus lead to vanishing 4NFs. While 
this is indeed the case for the leading 3N diagrams at 
NLO, it is explicitly shown in Ref. (186; 192) that many 
of the reducible-like diagrams in Fig. 14 do generate non- 
vanishing 4NF's which are not suppressed by inverse pow- 
ers of the nucleon mass. As a representative example, we 
give here the g-contribution which results entirely from 
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the first diagram in Fig. 14 (the second graph appears to 
be truly reducible and does not produce any contribution 
to the 4NF): 


yvi® — 29a HG 

ad (2F)® (@? + M2] [a + M2]? [@? + MZ] 
x (ri tate+t3— Ti Ty 72°74) di Gitte 

+71 X7T2°TA GT Ti2 X Ga 73 

+71 X73°T4%° U2 NT X V2: Fe (2.25) 
+71°T4 Gi2 X + 2 Gi2 X G4 73) + all perm., 
where G2 = Gi + @ = -—% — @ = —G@z is the momen- 
tum transfer between the nucleon pairs 12 and 34. The 


complete expression for the leading 4NF both in momen- 
tum and coordinate space can be found in Ref. (186). 
A rough estimation of the 4NF contributions to e.g. the 
a-particle binding energy is provided by the strength of 
the corresponding r-space potentials expressed in terms 
of dimensionless variables r;;M,. One then finds e.g. for 
the g-terms g§ M7(167F?)~-3 ~ 50 keV. This agrees 
qualitatively with a more accurate numerical estimation 
carried out in Ref. (195) which, however, still involved 
severe approximations to simplify the calculations. 

So far we only discussed isospin-invariant contribu- 
tions to the nuclear forces. It is well established that 
nuclear forces are charge dependent (for reviews see 
e.g. (196; 197)). For example, in the nucleon-nucleon 
199 channel one has for the scattering lengths a and the 
effective ranges r (after removing electromagnetic effects) 


acip = = (@nn + App) — Anp = 5.64 + 0.40 fm , 


NLR wele 


roi = = (Tan +Tpp) — Tnp = 0.03 + 0.06 fm (2.26) 


These numbers for charge independence breaking (CIB) 
are taken from the recent compilation of Machleidt (18). 
The charge independence breaking in the scattering 
lengths is large, of the order of 25%, since @np = 
(—23.74040.020) fm. Of course, it is magnified at thresh- 
old due to kinematic factors (as witnessed by the dis- 
appearance of the effect in the effective range). In ad- 
dition, there are charge symmetry breaking (CSB) ef 
fects leading to different values for the pp and nn phase 
shifts/threshold parameters, 


QCSB = App — Ann = 1.6+ 0.6 fm , 


TCSB = Tpp —Tnn =0.10+0.12 fm. (2.27) 
Combining these numbers gives as central values any = 
—18.9fm and ap, = —17.3fm. Notice that this value 
for Gyn is in agreement with the recent experimental 
determinations from the reaction m~d — nny, ann = 
—18.5+0.5fm (198), and the kinematically complete 
deuteron breakup reaction nd — nnp at Fjap = 13 MeV, 
Qnn = —18.7+0.6fm (199). However, another recent 
experiment also based on the deuteron breakup reaction 


at Ejap = 25.2 MeV yielded a considerably smaller value, 
Qnn = —16.3+ 0.4fm (200). For a review of indirect 
methods to measure the 1So nn scattering length and 
the current experimental status for this observable see 
(201). 

Within the Standard Model, isospin violation has its 
origin in the different masses of the up and down quarks 
and the electromagnetic interactions. Chiral EFT is well 
suited to explore the consequences of these two effects for 
low-energy dynamics of few- and many-nucleon systems. 
Consider first the strong isospin-violating effects. The 
QCD quark mass term can be written in the two-flavor 
case as 


1 
Linas = — 5 F(Mu + ma)(1 + er") a, (2.28) 
where the superscript of the Pauli isospin matrix denotes 
the corresponding Cartesian component and 
My — Ma 1 


nw, 2.29 
My + Ma 3 ( ) 


€ 


Here, the numerical estimation corresponds to the mod- 
ified MS subtraction scheme at a renormalization scale 
of 1 GeV (202). In Eq. (2.28), the isoscalar term breaks 
chiral but preserves isospin symmetry. It is responsible 
e.g. for the nonvanishing pion mass M? = (m, + mq)B, 
M2 = M? + O(m2 ,), and generates a string of chiral- 
symmetry-breaking terms in the effective hadronic La- 
grangian proportional to M?” with n = 1,2,.... The 
isovector term gives rise to the strong isospin breaking 
and leads to hadronic effective interactions «x («M?)”. 
Consequently, the typical size of the strong isospin vio- 
lation in hadronic observables is given by «M?/A? ~ 1% 
if one takes A = M, (this, however, does not apply 
e.g. to the pion masses). The leading and subleading 
strong isospin-violating contributions are already incor- 
porated in the Lagrangians £®) and gu in Eqs. (2.9), 
(2.12), (2.17) and correspond to terms involving y, +. 
Notice that the strong isospin violation is additionally 
suppressed in the meson sector (due to G-parity). In 
particular, the charged-to-neutral pion mass difference is 
almost entirely of electromagnetic origin. Electromag- 
netic terms in the effective Lagrangian resulting from ex- 
change of hard virtual photons can be generated using 
the method of external sources (35). All such terms are 
proportional to positive powers of the nucleon charge ma- 
trix Q = e(1+7°)/2 where e denotes the electric charge. 
In addition, soft photons have to be included explicitly. 
For more details on the inclusion of virtual photons in 
chiral EFT the reader is referred to (203-209). 

To explore isospin-breaking (IB) effects in nuclear 
forces/few-nucleon observables it is useful to relate the 
corresponding small parameters € and e to the chiral ex- 
pansion parameter M,/A. Clearly, this can be done in 
various ways. For example, in Refs. (210-212) the fol- 
lowing rules have been adopted: 

M, e? M4 
Ne np oe ee (2.30) 


Eere~ 


19 


Similar counting rules were also used in Refs. (213-218). 
Notice, however, that in the meson and single-baryon sec- 
tors one usually counts « ~ 1 but e~ M,/A. Utilizing 
the counting rules in Eq. (2.30), the leading and sublead- 
ing IB contributions from the hard virtual photons have 
the form (208): 


£2 em = C(QUQU), 
LEN em = FON(f(O2 - 2) + folQ4)O+ 
+ f3(Q2+Q2))N, (2.31) 
where Q4 = — 1/2(Qx) and f; refer to the corre- 


sponding LECs. The leading Lagrangians for the strong 
and electromagnetic IB contact interactions Les str ond 
LO). an are given explicitly in Ref. (210). 

Let us first discuss IB contributions to the pion and 
nucleon masses. As already pointed out, the leading con- 
tribution to the charged—to—neutral pion mass difference 
is entirely of electromagnetic origin, 


2 
6M2 = M2, —- M2, ~ eC. 


a (2.32) 


The experimentally known pion mass difference M,+ — 
M,0 = 4.6 MeV allows to fix the value of the LEC C, 
C = 5.9-10-° GeV*. Notice that the natural scale for 
this LEC is F?A?/(47)? ~ 3-107-° GeV“ if one adopts 
Aw~ M,. Writing the nucleon mass m as 


1 
m = diag(m,, mn) =m+ somr , (2.33) 


one obtains for the proton-to-neutron mass difference 6m: 


6m = —4cseM? — foe? F24+..., (2.34) 
where the ellipses refer to higher-order corrections. No- 
tice that the f3-term in Eq. (2.31) is isospin-invariant. 
On the other hand, the fi-term does produce IB ver- 
tices with two and more pions but does not contribute 
to the leading electromagnetic nucleon mass shift. The 
LECs cs and fg can be determined from the strong and 
electromagnetic nucleon mass shifts: 


(Mp — Mn)" = (dm)** = —2.05 + 0.3 MeV , 
(Mp — My)" = (6m)°*™ = 0.76 + 0.3 MeV , (2.35) 


which leads to (207) 


C5 = —0.09+ 0.01 GeV, 


fo = —0.45 £0.19 GeV? 


The values for the strong and electromagnetic nucleon 
mass shifts are taken from (219). Notice that the recent 
lattice QCD result (220) for the strong nucleon mass shift 
(dm)s* = —2.26 + 0.57 + 0.42 + 0.10 MeV is in agree- 
ment with the one of Ref. (219). We further emphasize 
that according to the counting rules in Eq. (2.30), the 
electromagnetic contribution to the nucleon mass shift is 
formally of higher order than the strong one. Based on 


naive dimensional analysis, these contributions are ex- 
pected to be of the size |(5m)*"| ~ eM?2/M, ~ 8 MeV 
and |(6m)*™| ~ e?M,/(47)? ~ 0.5 MeV. 

We are now in the position to overview the structure of 
the IB nuclear forces. The general isospin structure of the 
two-nucleon force feeds, according to the classification of 
Henley and Miller (221), into the four classes: Viy = 
a+ 71-72 (isospin-invariant), Vii = arr} (charge- 
independence-breaking), Vi! = a(7} + 73) (charge- 
symmetry-breaking) and Viv = 6 = a (73-73) +8 [71 x Te]? 
(issospin-mixing). Here, a anil GB denote the correspond- 
ing space and spin operators. Notice that for the class-IV 
terms, 3 has to be odd under a time reversal transforma- 
tion. The most general isospin structure of the 3NF is 
worked out in Ref. (211). While the distinction between 
the class-I, II and III forces based on the conservation 
of the total isospin operator T = ()),7;)/2 and charge- 
symmetry operator P., = exp(i7T2) can be straightfor- 
wardly generalized to any number of nucleons, the conser- 
vation of the operator T” responsible for the distinction 
between the class-III and IV 2NFs depends, in general, 
on the number of nucleons. In particular, the class-II 
and III 2NFs commute with the operator Ty (i.e. do 
not mix isospin in the 2N system) but do not commute 
with Te n- For this reason, the general isospin structure 
of the 3NF was classified in Ref. (211) in terms of class-I, 
II and III contributions. 

The dominant IB contribution to the 2NF occurs at 
NL@ (the slash indicates that we now use the power 
counting rules extended as in Eq. (2.30)) due to the 
charge-to-neutral pion mass difference in the 17-exchange 
potential. It can be accounted for by taking the proper 
pion masses in the la-exchange potential for various 
physical channels: 


vip = viet 
Ls es 
Vie _ 


= Vin(M, 0) ) 


—Vin(M. a) + 2(- 1)'t1Vi (M+) ’ (2.36) 
where J denotes the total isospin of the two-nucleon sys- 
tem and 
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(2.37) 


Notice that the resulting IB interaction conserves charge 
symmetry (i.e. class-II) and reaches about 6M?/M?2 ~ 
7% of the strength of the isospin-invariant 17-exchange 
potential. It is known to yield a sizeable contribution to 
the CIB in the 1So NN scattering length, see e.g. (210). 
Another IB effect at the same order comes from the 
Coulomb interaction between the protons (class-II and 
III). We emphasize that effects of the purely electromag- 
netic interactions in two-nucleon scattering observables 
get enhanced under certain kinematical conditions (low 
energies and/or forward angles) due to the long-range 
nature of these interactions. Clearly, such an enhance- 
ment goes beyond the simple power counting rules in 
Eq. (2.30). Consequently, despite the fact that the first 
corrections to the point-like static one-photon exchange 
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(Coulomb interaction) due to recoil and two-photon ex- 
change (222), pion loop contributions to the nucleon form 
factors (223), vacuum polarization (224; 225) and mag- 
netic moment interaction (226) are suppressed by the 
factor 1/m? and, according to the power counting, con- 
tribute at rather high orders, sizeable effects may show 
up in certain observables. For example, the magnetic 
moment interaction strongly affects the nucleon analyz- 
ing power A, at low energy and forward angles. Effects of 
subleading electromagnetic interactions were also inves- 
tigated in 3N continuum using phenomenological nuclear 
forces (227; 228), see (229) for a formulation based on 
pionless EFT. 


The corrections to the IB 2NF at N*L@ are CSB 
and arise from charge dependence of the pion-nucleon 
coupling constant in the 1la-exchange potential and 
the derivative-less NN contact interaction x my, — mq 
(212; 214). Notice, however, that the energy-dependent 
Nijmegen PWA does not yield any evidence for charge 
dependence of the pion-nucleon coupling constant (230). 
The leading CIB contact interactions are of electromag- 
netic origin and (formally) start to contribute at N°L@. 
At this order, one also has to take into account fur- 
ther IB contributions to the lz-exchange potential at 
the one-loop level which, to a large extent, can be ac- 
counted for by a further (charge-dependent) renormal- 
ization of the aN coupling constants in Eq. (2.36). The 
only contributions which have a different momentum de- 
pendence and, therefore, cannot be cast into the form 
of Eq. (2.36) are the ones « (6M,)? and the proton- 
to-neutron mass difference which involve class-IV oper- 
ators (212; 217), see (231) for a much earlier derivation 
of these terms. Notice that the power counting rules in 
Eq. (2. 30) suggests the following hierarchy of the 2NF 
(213): Viéy > Vii > VA > VEY which is consis- 
tent with the Seen acne: Next: my-exchange also con- 
tributes at this order. The resulting CIB potential has 
been worked out by van Kolck at al. (232) and re-derived 
recently by Kaiser (223). It can be written in a rather 
compact way and leads to negligibly small effects in NN 
scattering. Kaiser also calculated subleading contribu- 
tions to the my-exchange potential proportional to the 
large isovector magnetic moment k, = 4.7 of the nu- 
cleon and found that the resulting potentials, which are 
also CIB, have a similar strength as the leading-order 
one (223).#4 IB 27-exchange also starts to contribute at 
N°LQ and is driven by the neutral-to-charged pion mass 
difference (CIB) (215) and the strong contribution to the 
nucleon mass shift (CSB) (212; 216; 233; 234), see (210) 
for the application to NN phase shifts. Finally, there 
are also the first IB 3NFs. While the dominant CIB 27- 
exchange 2NF is generated by the pion mass difference, 


#4 Notice, however, that these corrections are suppressed by the 
factor 1/m relative to the leading-order contributions and, there- 
fore, appear formally at NOL@. 


the 3N diagrams with one insertion of 5M? at N3L@ are 
additionaly suppressed by the factor 1/m (if one uses 
an energy-independent formulation), see the discussion 
about the 3NFs at NLO earlier in the text. The nonva- 
nishing 3NF's at N°L@ result from 17- and 27-exchange 
3N diagrams constructed with the leading-order isospin- 
invariant vertices and a single insertion of 6m as well 
as 27-exchange diagram with the leading IB m7 NN in- 
teractions « fi,2 (211; 218). One finds that all these 
contributions are CSB except the one which is propor- 
tional to LEC f; and is CIB. We further emphasize that 
while the value of the LEC f2 is determined by the elec- 
tromagnetic nucleon mass shift, the value of the LEC fy; 
is unknown. However, see (209) for an estimation of fi 
based on dimensional analysis and (235) for an attempt 
to determine f; from data. 


Remarkably, even the N*LQ contributions to the two- 
and three-nucleon forces have been worked out. At this 
order, no new structures appear in the 17-exchange po- 
tentials. The corrections to the leading IB 27-exchange 
potential result from a single insertion of either the sub- 
leading isospin-conserving m7NWN vertices proportional 
to the LECs c;, see Eq. (2.17), the leading electromag- 
netic vertex proportional to the LEC fe, see Eq. (2.31), 
#5 or the (poorly known) leading charge-dependance of 
the pion-nucleon coupling constant (212). The resulting 
IB potentials involve the class-II and III central, tensor 
and spin-spin components. The CIB potentials typically 
have the strength of a few 10 keV at relative distances 
r ~ M;'. The CSB tensor and spin-spin potentials are 
weaker (< 10 keV) while the CSB central potential is 
comparable in size to the CIB contributions. Similarly 
to the isospin-conserving 27-exchange potential, the sub- 
leading contributions turn out to be numerically large in 
comparison to the leading-order ones. In particular, for 
the class-III central 27-exchange potential one obtains 
eon () ~ 3 for r ~ Mz}. The main reason 
for this unpleasant convergence pattern is the same as 
in the isospin-conserving case and can be traced back to 
the (large) A-isobar contributions to the LECs c3,4. We 
will discuss this issue in more detail in section II.D. Last 
but not least, there are also numerous IB contact inter- 
actions with up to two derivatives involving class-II, III 
and IV terms, see also (217). The corrections to the 3NF 
at N4L@ are worked out in Refs. (211; 218). At this or- 
der, the first IB but charge symmerty conserving 3NFs 
show up which result from the neutral-to-charged pion 
mass difference in the N?LO 2z7- and 17-exchange dia- 
grams in the second line of Fig. 13 and the 27-exchange 
diagram involving the m7NN vertex « fi. In addi- 
tion, there are CSB 3NFs of the 27- and 17-exchange 
types driven by the electromagnetic nucleon mass shift. 
Again, the strongest 3NF's turn out to be the ones which 


#5 The two other LECs f1,3 do not contribute to the 27-exchange 
2NF at this order. 
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are proportional to the LECs c34. They are charge- 
symmetry conserving and arise from a single insertion 
of 6M? into the pion propagators of the 27-exchange 3N 
graph. The expected strength of such IB potentials is 
~ 26M?/M?2 ~ 13% as compared to the isospin-invariant 
ones given in Eq. (2.22) which are known to yield about 
~ 500...1000 keV to the triton binding energy (the 
precise numbers are renormalization scheme dependent). 
Also the strength of the corresponding coordinate-space 
potentials, e.g. |g36M?2M4c3/(641?F+)| ~ 70 keV (here 
we picked out one particular term), indicates that the 
resulting IB effects in few-nucleon observables might be 
sizeable. The CSB 3NFs, on the other hand, do show a 
more natural convergence pattern and are considerably 
waeker. Their contribution to e.g. the 3H-?He binding 
energy difference is expected to be of the order ~ 10 
keV. 

Recently, a certain classes of even higher-order con- 
tributions have been worked out by Kaiser. In particu- 
lar, he calculated the subleading w77y-exchange potentials 
proportional to the LECs c34 at the 2-loop level which 
(formally) contribute at order N°L@ (236; 237). Espe- 
cially the contributions driven by the LEC c3 were found 
to generate astonishingly strong CSB and CIB potentials 
which amount to ~ 1% of the strongly attractive isoscalar 
central potential at N?LO and reach a few 100s of keV 
at r ~ M~'. Notice, however, that effects of these very 
strong potentials in S-, P- and D-waves may, to some 
extent, be compensated by the corresponding IB con- 
tact interactions. The effects in higher partial waves are 
presumably suppressed due to the shorter range of the 
27m-exchange potential compared with the 17-exchange 
one. 


C. Chiral EFT for few nucleons: applications 


We now turn our attention to applications. As dis- 
cussed in the previous section, the two-nucleon chi- 
ral potential involves the long-range contributions due 
to the multiple pion exchanges and short-range ones 
parametrized by contact interactions. Both kinds of 
terms typically grow with increasing nucleon momenta 
and become meaningless in the large-momentum re- 
gion as follows from the very nature of EFT being 
the low-momentum expansion. As a consequence, the 
Schrodinger equation is ultraviolet divergent and needs 
to be regularized (and renormalized). The problem of 
renormalization in the nonperturbative regime in the con- 
text of both pion-less (41; 238-245) and pion-full EFT 
(139; 140; 168; 246-266) has attracted a lot of interest in 
the past years. The standard procedure to renormalize 
the Lippmann-Schwinger (LS) equation is based on Wil- 
son’s method and implies the following two steps (246). 
First, one solves the LS equation regularized with the 
finite momentum (or coordinate-space) cutoff and using 
the potential truncated at a given order in the chiral ex- 
pansion as the kernel. Secondly, the LECs accompany- 


ing the contact terms in the potential are determined 
by matching the resulting phase shifts to experimental 
data which, in this framework, can be viewed as renor- 
malization. Notice that iterating the truncated expres- 
sion for the chiral potential in the LS equation neces- 
sarily generates ultraviolet divergencies in the Neumann 
series which require counterterms beyond the given ap- 
proximation for the potential. As a consequence, taking 
the limit of the infinite cutoff in such a manifestly non- 
renormalizable (in the above mentioned sense) approach 
might result e.g. in impossibility to resolve the (nonlin- 
ear) matching conditions for the corresponding LECs. 
A detailed discussion on the choice of ultraviolet cutoff 
and its role in renormalization of the Schrodinger equa- 
tion is given by Lepage (246; 247). He argued that the 
coordinate-space (momentum-space) cutoff should not be 
decreased (increased) beyond the separation scale, after 
which the description of the data stops to improve. Tak- 
ing the cutoff near this separation scale is the most ef- 
ficient choice. This strategy has been followed by the 
currently most advanced N?LO analyses of the 2N sys- 
tem of Refs. (271; 272) where the cutoffs A = 450...600 
MeV have been employed. These studies were critisized 
by Nogga et al. (139) who considered low NN partial 
waves based on the la-exchange potential and contact 
interactions using a much bigger cutoff variation with 
A <4 GeV. They found that higher-order counterterms 
have to be promoted to LO in the °Po, ?Ps-°F> and 2D» 
channels in order to stabilize the amplitude. On the other 
hand, the efficiency of such a modified power counting 
framework was questioned in Ref. (258), where it has 
been demonstrated that increasing the cutoff and pro- 
moting counterterms as suggested in Ref. (139) does not 
improve the overall description of the scattering observ- 
ables. For more discussions on the conceptual issues re- 
lated to the power counting in the NN system the reader 
is referred to (139; 140; 168; 246; 247; 254; 258; 262). 
More work is needed in the future in order to clarify the 
relation between the well-established chiral expansion of 
the nuclear potential and the scattering amplitude. 


We further emphasize that it is possible to nonpertur- 
batively renormalize the partial-wave-projected LS equa- 
tion with singular 1/r”-potentials (112; 262; 267-270). 
This program was applied to different NN channels based 
on the la- and 27-exchange potentials at various or- 
ders in the chiral expansion by the Granada group 
(252; 253; 255-257; 260; 263; 264). In these studies, 
the short-range counter terms are replaced by adjustable 
parameters entering the short-distance boundary condi- 
tions. The number of such parameters in each channel is 
uniquely determined by the sign (attractive vs. repulsive) 
of the strongest singularity which raises concerns about 
a systematic improvability (in the EFT sense) of such a 
framework. Nevertheless, the findings of these studies in 
attractive channels provide an impressive demonstration 
of the existence of the long-range correlations in the NN 
scattering observables. 


The most advanced analyses of the two-nucleon sys- 
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FIG. 15 Neutron-proton phase shifts in S- and P-waves at 
N8LO in comparison with the Nijmegen (273; 274) (filled cir- 
cles) and Virginia Tech (275) (open triangles) PWA. Shaded 
bands (dashed lines) refer to the calculations by EGM (272) 
(EM (271)). 


tem based on the Weinberg power counting take into 
account the 2NF contributions up to N°LO (271; 272). 
Most of the LECs c;, dj entering the long-range part of 
the potential are sufficiently well determined in the pion- 
nucleon system (184).#° The 24#7 unknown LECs en- 
tering the short-range part of the 2NF at N°LO have 
been extracted from the low-energy NN data for several 
choices of the cutoff in the Schrédinger equation. Both 
N°LO potentials of Entem and Machleidt (271) (EM) 
and Epelbaum, Gloéckle and MeifSner (272) (EGM) yield 
accurate results for the neutron-proton phase shifts up 
to Flap ~ 200 MeV and the deuteron observables. This 
is exemplified in Figs. 15, 16 where the EGM and EM 
results for the neutron-proton S- P- and D-waves and 
the corresponding mixing angles are shown in compari- 
son with PWA results from Refs. (273-275). The bands 
in the EGM analysis result from the variation of the cut- 
off in the LS equation (spectral function regularization) 
in the range A = 450...600 MeV (A = 500...700 MeV). 
It is comforting to see that in most cases the results of 
both analyses agree with each other within the estimated 
theoretical uncertainty. Notice, however, that the EM 


#6 Notice, however, that the value of the LEC ca adopted in (271), 
c4 = 5.4 GeV—! is not compatible with pion-nucleon scattering 
where one finds at order Q? (185): cg = 3.40 £0.04 GeV-!. 

#7 This number refers to isospin-invariant contact interactions. 
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FIG. 16 Neutron-proton phase shifts in D-waves and the mix- 
ing angles €1,2 at N°LO. For notation see Fig. 15. 


and EGM analyses differ from each other in several im- 
portant aspects. For example, the so-called spectral func- 
tion regularization (155; 156) of the 27-exchange contri- 
butions has been adopted by EGM while the analysis by 
EM is based upon dimensionally regularized expressions. 
Further differences can be attributed to the implemen- 
tation of the momentum-space cutoff in the Schrédinger 
equation and the treatment of relativistic effects. More 
precisely, the work by EGM is based on the “relativis- 
tic” Schrodinger /Lippmann-Schwinger equation, a natu- 
ral extension of the usual nonrelativistic equations uti- 
lizing the relativistic relation between the CMS energy 
and momentum, see (170) for more details. This equa- 
tion can be straightforwardly generalized to the case of 
several nucleons, see (276-278) for recent studies of rela- 
tivistic effects in 3N observables, and can also be cast into 
equivalent nonrelativistic-like forms (170; 279) (provided 
the potential is appropriately modified). On the other 
hand, the analysis of Ref. (271) is based on the non- 
relativistic Schrodinger equation and uses the static I7- 
exchange potential and the 1/m- and 1/m?-corrections to 
the 27-exchange potential from Refs. (142; 173), where 
no particular dynamical equation is specified. It is un- 
clear whether the relativistic corrections to the poten- 
tial used in (271) are consistent with the nonrelativistic 
Schrodinger equation. Further differences between the 
EGM and EM analyses result from the fitting procedure: 
the LECs accompanying the contact interactions were 
determined by EM/EGM by fitting directly to the scat- 
tering data/to the Nijmegen PWA (273; 274). For this 
reason, EGM adopted the same treatment of IB effects 
as followed by the Nijmegen group and did not include 
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FIG. 17 Neutron-proton So partial wave at NLO (dashed 
band), N?LO (light-shaded band) and N?LO (dark-shaded 
band) in comparison with the Nijmegen (273; 274) (filled cir- 
cles) and Virginia Tech (275) (open triangles) PWA. 


e.g. the leading IB contributions to the 27-exchange po- 
tential. Perhaps, the most important difference between 
the two studies is related to the estimation of the theo- 
retical uncertainty. In the work by EGM, the theoreti- 
cal uncertainty was estimated by varying the cutoffs in 
the Schrédinger equation and the spectral function repre- 
sentation for the 27-exchange potential restricted by the 
condition that the resulting LECs are of a natural size 
which might be viewed as a self-consistency check for cal- 
culations carried out within the power counting scheme 
based on naive dimensional analysis, see (87; 163). No- 
tice further that at NLO and N?LO, the strengths of 
various contact interactions are well understood in terms 
of resonance saturation on the basis of phenomenological 
one-boson exchange models (280). No serious attempt 
to provide a realistic error estimation was done in the 
analysis of EM. On the other hand, their work clearly 
demonstrates that for a particularly chosen regulariza- 
tion prescription it is even possible to accuratly describe 
two-nucleon scattering data for Ejay > 200 MeV. For 
further technical details, results for various scattering ob- 
servables and the properties of the deuteron the reader 
is referred to the original publications (271; 272) and the 
review article (163). 


To illustrate the convergence of the chiral expansion 
for NN phase shifts, we show in Fig. 17 the results for 
the So partial wave at NLO, N?LO (156) and N3LO 
(272). We emphasize that the variation of the cutoff at 
both NLO and N?LO only shows the effects of missing 
N°LO contact interactions. It, therefore, does not pro- 
vide a realistic estimation of the theoretical uncertainty 
at NLO, see (163) for an extended discussion. 


Applications to the three-nucleon system have so far 
been carried out up to N2LO. At NLO, no 3NF needs 
to be taken into account. This allowed for a parameter- 
free predictions of various 3N scattering observables at 
low energies as well as for the triton and a-particle bind- 
ing energies (122). Using the most recent version of the 


NLO potential based on the spectral functions regular- 
ization, one finds at NLO (163) Bsy = 7.71...8.46 MeV 
and Bsy. = 24.38...28.77 MeV to be compared with 
the experimental values Bsy = 8.482 MeV and Bay, = 
28.30 MeV. These numbers are similar to the ones ob- 
tained in Ref. (122) within the framework based on di- 
mensional regularization. 

At N?LO one, for the first time, has to take into ac- 
count the crresponding 3NFs. The two LECs D and E 
entering the expressions for the 3NF in Eq. (2.22) have 
been determined by fitting the °H binding energy and 
either the nd doublet scattering length (123), the *He 
binding energy (281) or the properties of light nuclei 
(282). Notice that the 7NNNWN vertex entering the 
17-exchange-contact 3NF also plays an important role 
in processes with a completely different kinematics such 
as e.g. the pion production in the NN collisions (283), 
see section II.E, or weak reactions like pp > deTv., see 
(284) and references therein. This offers the possibil- 
ity to extract the corresponding LEC from these pro- 
cesses, see (284) for a recent attempt. With the LECs 
being determined as described above, the resulting nu- 
clear Hamiltonian can be used to describe the dynam- 
ics of few-nucleon systems. In particular, 3N continuum 
observables offer a natural and rich testing ground for 
the chiral forces. In Refs. (122; 123; 285-293) various 
3N scattering observables have been explored by solving 
the momentum-space Faddeev equations with chiral two- 
and three-nuclein forces as input. In the formulation of 
Ref. (22), one first computes the T-matrix by solving the 
Faddeev-like integral equation 


T = tPd+(1+tGo) Vey (1+ P)¢+tPGoT 
+ (1+tGo) Vey 1+ P)GoT, (2.38) 


where the initial state ¢@ is composed of a deuteron and 
a momentum eigenstate of the projectile nucleon. Here 
V3 is that part of the 3N force which singles out the par- 
ticle 2 and which is symmetric under the interchange of 
the two other particles. The complete 3NF is decomposed 
as Van = V3iy + Vey + Vey. Further, Go = 1/(E — Ho) 
is the free propagator of the nucleons, P is a sum of a 
cyclical and anti-cyclical permutation of the three par- 
ticles and t denotes the two-body t-matrix. Once T is 
calculated, the transition operators U.; and U,, for the 
elastic and break-up channels can be obtained via 


Ua = PGo'+PT+Vay(1+P)(1+GoT), 
Up = (1+P)T. (2.39) 


For details on solving these equations in momentum 
space using a partial wave decomposition the reader is 
referred to (294). The partial wave decomposition of 
the 17-exchange and contact 3NF at N?LO and the one- 
pion-two-pion-exchange topology at N°LO is detailed in 
Refs. (123) and (295), respectively. The expressions for 
various observables in terms of the transition operators 
are given in (22). The inclusion of the long-range electro- 
magnetic interaction requires a non-trivial generalization 
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FIG. 18 Differential cross section for elastic nd scattering at 
E\ab = 10 MeV (left panel) and 65 MeV (right panel). Light 
(dark) shaded bands depict the results at NLO (N?LO). The 
neutron-deuteron data at 10 MeV are from Ref. (298). The 
remaining data at 10 MeV are the Coulomb/IB-corrected 
proton-deuteron data from Refs. (299-301). The data at 
65 MeV are proton-deuteron data from Ref. (302). 


of the formalism, see (296; 297) for recent progress along 
this line. 


The results for the differential cross section in elas- 
tic nd scattering are in a good agreement with the data, 
see Fig. 18 for two representative examples. Notice, how- 
ever, that the theoretical uncertainty becomes significant 
already at intermediate energies. Qualitatively, this be- 
havior is consistent with the one observed in the two- 
nucleon system (272). Notice further that the descrip- 
tion of the data improves significantly when going from 
NLO to N?LO. The situation is similar for vector and 
tensor analyzing powers, see Ref. (163) for a recent re- 
view article. More complicated spin observables have also 
been studied. As a representative example, we show in 
Fig. 19 a selection of the proton-to-proton and proton- 
to-deuteron polarization transfer coefficients measured in 
d(p, p)d and d(p, d)p reactions at EG = 22.7 MeV 
(303; 304). The results at N?LO are in a reasonable 
agreement with the data, see (302) for more examples. 
One further observes that the theoretical uncertainty ob- 
tained by the cutoff variation is underestimated at NLO, 
see the discussion earlier in the text. It is, however, com- 
forting to see that the description of the data improves 
significantly when going from NLO to N?LO. 


The nucleon-deuteron breakup reaction offers even 
more possibilities than the elastic channel due to the 
much richer kinematics corresponding to three nucleons 
in the final state. It has also been studied extensively 
over the last years, both theoretically and experimentally, 
leaving one with mixed conclusions. While the differen- 
tial cross section in some configurations such as e.g. the 
recently measured np final-state interaction, co-planar 
star and an intermediate-star geometries at low energies 
are in a very good agreement with the data (287), large 
deviations are observed in certain other configurations. 
In particular, the so-called symmetric space-star configu- 
ration (SST) appears rather puzzling. In this configura- 
tion, the plane in the CMS spanned by the outgoing nu- 
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FIG. 19 The proton-to-proton (left panel) and proton-to- 
deuteron (right panel) polarization transfer coefficients in 
d(p, ~)d and d(p, d)p reactions at Ee? = 22.7. Light (dark) 
shaded bands depict the results at NLO (N?LO). Data are 
from Refs. (303; 304). 
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FIG. 20 Chiral EFT predictions for neutron-deuteron 
breakup cross section (in mb MeV~! sr~*) along the kinemat- 
ical locus S. Light-shaded (dark-shaded) bands refer to the 
results at NLO (N?LO). Left panel: The SST configuration 
at En = 13 MeV. Neutron-deuteron data (open triangles) are 
from (305; 306), proton-deuteron data (filled circles) are from 
(307). Right panel: The SCRE configuration with a = 56° 
at EN = 19 MeV (292). Dashed and dashed-dotted lines are 
results based on the CD Bonn 2000 2NF (18) combined with 
the TM99 3NF (308) and the coupled channel calculation in- 
cluding the explicit A and the Coulomb interaction (296), 
respectively. 


cleons is perpendicular to the beam axis, and the angles 
between the nucleons are 120°. At Ejay = 13 MeV, the 
proton-deuteron and neutron-deuteron (nd) cross section 
data deviate significantly from each other. Theoretical 
calculations based on both phenomenological and chiral 
nuclear forces have been carried out for the nd case and 
are unable to describe the data, see Fig. 20. Moreover, 
the Coulomb effect was found to be far too small to ex- 
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plain the difference between the pd and nd data sets. 
Recently, proton-deuteron data for a similar symmetric 
constant relative-energy (SCRE) configuration have been 
measured in Cologne (292). This geometry is character- 
ized by the angle a between the beam axis and the plane 
in the CMS spanned by the outgoing nucleons. Similar 
to the SST geometry, one observes large deviations be- 
tween the theory and the data, in particular for a = 56°, 
see Fig. 20. The included 3NFs have little effect on the 
cross section while the effect of the Coulomb interaction 
is significant and removes a part of the discrepancy. No- 
tice that all above cases correspond to rather low energies 
where one expects good convergence of the chiral expan- 
sion. Furthermore, contrary to the A,-puzzle, the cross 
sections discussed above are mainly sensitive to the two- 
nucleon S-waves without any known fine tuning between 
partial waves. First attempts have been made in the 
past few years to perform deuteron breakup experiments 
at intermediate energies, in particular at Ey = 65 MeV 
(289), in which a large part of the phase space is covered 
at once. Chiral EFT results at N?LO for more than 155 
data points were shown to be of a comparable quality 
to the ones based on modern phenomenological nuclear 
forces. 

Recently, first results for the 4N continuum based on 
both phenomenological and chiral nuclear forces and in- 
cluding the Coulomb interactions have become available, 
see (309; 310) for p—*He scattering, (311) for the n—*He, 
p—H and d—d scattering, and (312) for the related ear- 
lier work. These studies do not yet include effects of 
3NFs but clearly indicate that at least some of the puz- 
zles observed in the 3N continuum also persist in the 4N 
continuum (such as e.g. the A,-puzzle in p—*He scatter- 
ing (310)). For a promising new approach to describe 
scattering states in even heavier systems the reader is 
referred to (313). 

The properties of certain S-shell and P-shell nuclei 
with A < 13 have been analyzed recently based on the 
no-core shell model (NCSM), see (281; 282) and (314) 
for an overview. In Fig. 21 we show some results from 
Ref. (282) for the spectra of !°B, 11B, *C and 18C. We 
emphasize that the LECs D and E entering the N?LO 
3NF were determined in these calculations by the triton 
binding energy and a global fit to selected properties of 
STi, }°B and !2C. These studies clearly demonstrate that 
the chiral 3NF plays an important role in the descrip- 
tion of spectra and other properties of light nuclei. The 
inclusion of the 3NF allows to considerably improve the 
agreement with the data. Further results for light nu- 
clei and the dilute neutron matter based on the lattice 
formulation of chiral EFT are given in sections II.G and 
IILE. 


D. The role of the A-isobar 


The chiral expansion for the long-range part of the 
nuclear force discussed in the previous section exhibits a 
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FIG. 21 States dominated by P-shell configurations for ‘°B, 
"B,C and '3C. The excitation energy scales are in MeV. 
Calculation is carried out in the framework of NCSM based 
on chiral N°LO 2NF of Ref. (172) and N?LO 3NF. For more 
details on the calculation see (282). Figure courtesy of Petr 
Navratil. 


somewhat unnatural convergence pattern in certain cases 
such as e. g. for the central part of the 27-exchange po- 
tential. The origin of the unnaturally strong subleading 
contribution in this case potential can be traced back 
to the large values of the dimension-two low-energy con- 
stants (LECs) c3,4 which are also responsible for the nu- 
merical dominance of the subleading 37-exchange (180) 
and charge-symmetry breaking 27-exchange 2NF (212) 
over the corresponding leading contributions. The large 
values of these LECs are well understood in terms of res- 
onance saturation (160). In particular, the A(1232) pro- 
vides the dominant (significant) contribution to c3 (c4). 
Given its low excitation energy, A = ma—m = 293 MeV, 
and strong coupling to the 7N system, the A-isobar is 
known to play an important role in nuclear physics. One 
can, therefore, expect that the explicit inclusion of A in 
EFT will allow to resum a certain class of important con- 
tributions and improve the convergence as compared to 
the delta-less theory, provided a proper power counting 
scheme such as the small scale expansion (SSE) (315) is 
employed, see also (316). The SSE is a phenomenologi- 
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cal extension of chiral perturbation theory in which the 
small expansion parameter includes external momenta, 
pion masses and the nucleon-delta mass splitting, 


Q/A € {p/A, M/A, A/A} , 


i.e. the delta-nucleon mass splitting is treated as A = 
ma—m~ M, rather than A ~ A >> M,. The improved 
convergence has been explicitly demonstrated for pion- 
nucleon scattering where the description of the phase 
shifts at third order in the SSE comes out superior (in- 
ferior) to the third (fourth) order chiral expansion in the 
delta-less theory (319). In the following, we will overview 
the additional contributions to the nuclear force which 
arise in the A-full theory as compared to the A-less the- 
ory and discuss the implications for the convergence of 
the low-momentum expansion. Notice that in such a 
setting we do not need to consider the NN — NA, 
NN — AA and NA => AA transitions which would cor- 
respond to the coupled-channel approach. The adopted 
counting rules for the nucleon mass and the delta-nucleon 
mass difference imply for the momentum scale associated 
with real delta production ~mA ~ /mM, ~ A. For 
typical external momenta (energies) of the nucleons of 
the order |p| ~ M, (Exin ~ M?/m) we are interested 
in, the momenta associated with real delta production 
can be safely integrated out. We, therefore, only need 
to consider the contributions to the nuclear force arising 
due to virtual delta excitations. 

The effective Lagrangian can be straightforwardly ex- 
tended to include the A-degrees of freedom. To work 
out A-contributions up to N?LO, the following additional 
terms in the heavy-baryon Lagrangian have to be taken 
into account: 


(2.40) 


£2 = —Tiliv- DY — As] g’ Ti +..., 
Pa he TPH uN + he., 
Lora = (bs + bg) TeiPP’ wh}, v°N + hic. +.. (3.41) 


where T with y (i) being the Lorentz (isospin) index 
denotes the large component of the delta field. Fur- 
ther, D“! refers to the chiral covariant derivative for the 
delta fields and P,,, is the standard projector on the 3 /2- 
components, Puy = Juv — Upvy — 45,S,/(1 — d), with 
d being the number of space-time dimensions. We also 
have wi = (T'uq)/2 and wig = (7*[Aa, ua])/2. The only 
relevant LEC in the lowest-order Lagrangian is the 7NA 
axial coupling h4. At subleading order, the combination 
of zNA LECs b3 + bg contributes. For more details on 
the notation, the reader is referred to (315; 319), see also 
(38) for a recent review article and (317; 318) for differ- 
ent formulations. Finally, it should also be emphasized 
that the only possible derivative-less NN NA contact in- 
teraction 


LO”. « (THNNS,7'N +h. c.), (2.42) 


vanishes due to the Pauli principle (320). 


The values of the LECs in the 7N Lagrangian are, 
clearly, different in the A-less and A-full theories and 
can be naturally extracted from 7N scattering, see (319) 
for such a determination at the leading one-loop level 
(i.e. order Q*). At subleading order, which is sufficient 
for our purpose, the determination of c; from the 7N S- 
and P-wave threshold coefficients yields in the delta-less 
theory (321) 

cy = —0.57, co = 2.84, cz = —3.87, cq = 2.89, (2.43) 
where only central values are given and the units are 
GeV~-!. The above values are somewhat smaller in 
magnitude than the ones obtained at higher orders, see 


e.g. (184). Including the contributions form the A, one 
finds 


¢ = —0.57, c = —0.25, cs = —0.79, ca = 1.33, 
bebe =iA0: (2.44) 


Notice that the LECs c2.3,4 are strongly reduced in mag- 
nitude when the A-isobar is included. It should also be 
emphasized that the values of these LECs depend sen- 
sitively on the choice of h4, which in the above case 
was set to ha = 3g4/(2V2) from SU(4) (or large N-). 
The results for the threshold coefficients and the 27- 
exchange potential are, however, rather stable (321). We 
also emphasize that the description of the P-wave thresh- 
old parameters improves significantly upon inclusion of 
the delta-isobar. 

We are now in the position to discuss the leading and 
subleading contributions of the A-isobar to the nuclear 
force. Since the appearance of a virtual A-isobar requires 
at least one loop, the corresponding contributions first 
appear at NLO (v = 2). The relevant 2N and 3N di- 
agrams can be obtained from the ones of Figs. 12 and 
13 by replacing the nucleon propagators by the ones 
of the A-fields in all intermediate states. We first dis- 
cuss the 2NF. Similarly to the A-less theory, the ad- 
ditional contributions to the 1a-exchange potential and 
contact interactions at both NLO and N?LO only lead 
to renormalization of various LECs. The 27-exchange 
diagrams were first discussed by Ordénez et al. (150) 
using time-ordered perturbation theory. These contri- 
butions were then calculated by Kaiser et al. (322) us- 
ing the Feynman graph technique. The corrections at 
N?LO have been worked out recently (321). We refrain 
from showing here the resulting expressions which are 
rather involved and only give the results for the isovec- 
tor tensor 27-exchange potential W7, defined according 
to Van = 71-7201 -Go2-G Wr, which may serve as a 
representative example: 


h2, 

= A 
1296072 F4A 
+ (40? + w?)D¥(q))] , 


Wr = [9rg%w? AR (g) +n’ (2LA(q) 


27 


Wine [(2(bs + bs)ga(w? — 124%) 
TO  648r2F4 ls a 
— 9e4(w? — 4A?)) DA(q) 
+ 6(3e4 — 2(bs + bs)ha) L4(q)] . (2.45) 
Here, the new loop function DA(q) is defined via 
5 a a Jie — 42 
DAG Seo) Ls area = (AG) 
A Jom, MH’ +4 2A 


The complete results for the A-contributions can be 
found in Refs.(321; 322). It is instructive to verify the 
consistency between the A-full and A-less theories which 
requires that the contributions due to intermediate A- 
excitations, expanded in powers of 1/A can be absorbed 
into a redefinition of the LECs in the A-less theory. This 
is only possible if the nonpolynomial (in momenta) terms 
up to N?LO resulting from such an expansion have the 
same form as expressions in Eqs. (2.14, 2.18). This in- 
deed turns out to be the case: all expanded nonpolyno- 
mial terms up to N*LO are exactly reproduced by the 
shift in the LECs c3,4 


4h, 


ox (2.47) 


6 —2c4 = 
in Eqs. (2.14, 2.18). 

To get more insight into the strength of various 27- 
exchange contributions in the A-full and A-less theo- 
ries, it is useful to switch to coordinate space. The 27- 
exchange potential can then be written as 


V(r) i Vot71-T2Wet [Vs t71-T2 Ws] G1 + Ge 
+ [Vr +71 ‘TQ Wr] Sie, (2.48) 


where Siz = 301, -f O2 +f — 01 + G2 is the tensor op- 
erator. The scalar functions V;(r) and W;,(r) are plot- 
ted in Fig. 22 using the values for the LECs specified 
in Eqs. (2.43) and (2.44). As expected, one observes a 
more natural convergence pattern in the theory with ex- 
plicit deltas with the N?LO contributions yielding typ- 
ically only modest corrections to the NLO result. This 
is, clearly, not the case in the delta-less theory where 
the entire contributions to Vo and Wr. are generated 
at N?LO. On the other hand, the N?LO 27-exchange po- 
tential in the delta-less theory provides a very good ap- 
proximation to the potential resulting at the same order 
in the delta-full theory. This indicates that the satura- 
tion of the LECs c3,4 is the most important effect of the 
A-isobar at the considered order. The results for NN F- 
and other peripheral waves calculated using the Born ap- 
proximation also clearly demonstrate the improved con- 
vergence in the theory with explicit A, see Fig. 23. 

As explained in section II.B, the first nonvanishing 
contributions to the 3NF appear in the A-less theory 
at N?LO. The situation is different in the A-full theory 
where the first 3NF contribution is generated at NLO 
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FIG. 22 Isoscalar (left panel) and isovector (right panel) com- 
ponents of the 27-exchange potential in coordinate space for 


A = 700 MeV. Dashed and solid (dotted and dashed-dotted) 
lines refer to the NLO and N?LO results in the delta-full 
(delta-less) theory, respectively. There are no contributions to 
Vo and Wr,s (Vr,s and Wc) at NLO (N?LO) in the delta-less 
theory. 


by the second graph in the first line of Fig. 13 with the 
intermediate nucleon propagator being replaced by the 
one of the A-field. In fact, the importance of the 7N P33 
partial wave corresponding to the excitation of the A res- 
onance in the 3NF has been realized already 50 years by 
Fujita and Miyazawa (23). The resulting expression for 
the A-contribution to the 27-exchange 3NF is exactly re- 
produced by the first term in Eq. (2.22) if one uses the A- 
saturation values for the LECs c; from Eq. (2.47). This, 
in fact, follows from the decoupling theorem and the fact 
that the static A propagator is proportional to A~!. No- 
tice that there are no short-range 3NFs with intermediate 
A-excitation since the corresponding N N NA interaction 
is Pauli forbidden. Stated differently, the LECs D and 
F are not saturated by the A-isobar. Surprisingly, one 
finds that there are also no A-contributions to the 3NF 
at N?LO (320). The 27-exchange diagrams with one in- 
sertion of the subleading NA vertex « b3 + bg generate 
1/m-suppressed terms due to the time derivative enter- 
ing this vertex. Despite the fact that both the A-full and 
A-less theories yield the same expressions for the 27- 
exchange 3NF at N?LO, one should keep in mind that 
the strengths of various terms are different. The extrap- 
olation of the 7N amplitude from threshold, where the 
LECs are determined, to the kinematical region relevant 
for the 3NF is discussed by Pandharipande et al. (323) 
who claimed that both theories might yield sizeably dif- 
ferent results if the expansion is truncated at low orders. 
Using the values of the LECs from Eqs. (2.43) and (2.44) 
one, however, finds that the strengths of various terms 
differ from each other at most by ~ 7% at N?LO (320). 
To conclude, the only effect of including the A-isobar as 


FIG. 23 F-wave NN phase shifts for A = 700 MeV. The dot- 
ted curve is the LO prediction, long-dashed (short-dashed) 
and solid (dashed-dotted) lines show the NLO and N?LO re- 
sults with (without) the explicit A-contributions. The filled 
circles depict the results from the Nijmegen PWA (273). 


an explicit degree of freedom in the 3NF up to N?LO is 
the shift of the majour part of the 27-exchange contribu- 
tion in Eq. (2.22) from N?LO to NLO and some minor 
changes in the strengths of various terms in this expres- 
sion. 

We now discuss the role of the A for IB nuclear forces. 
The observed unnatural convergence pattern for the CSB 
27-exchange 2NF and CIB 3NF in the A-less theory, see 
the discussion in section II.B, is very similar to the one 
for isospin-conserving 27-exchange 2NF (in all cases the 
large contributions are proportional to the LECs cz,4) 
and provides a strong motivation to explore the role of 
the A isobar in this case. 

The leading IB A-contributions to the 27-exchange 
2NF result from the corresponding triangle, box and 
crossed-box diagrams with one insertion of isospin- 
breaking pion, nucleon and delta mass shifts. The latter 
can be deduced from the corresponding leading strong 
and electromagnetic Lagrangians (320) 


CoA = ~The (x4 — (x+)) 69 Suv TF 
Chm = TF EB big (Q4. — Q2) + f2 5i5 (Q4)Q4 
+ £3 biz (Q4)? + £2 (7°Q4)(77Q4) 
+ fS (r'Q_NriQ)] Gap TY Ban, (2:49) 
where cS and f* are the LECs and the ellipses in the last 
line refer to strong terms which involve at least one pion 


field and are irrelevant for the following discussion. The 
masses of the physical delta fields (Att, At, A°, A~) 


can be written as 


7 dmh z dm dm? 
Matt =Ma+— 4, mat =mat—44+—4, 
2 6 2 
_ bmi dm? 2 dm1 
mao = MA — 7 a MAaA- = Ma — 5 


where dm) /dm denote the equidistant /non-equidistant 
splittings and the mass m, contains an isospin-invariant 
shift dma defined as MA =m, +6ma, with ma being 
the delta mass in the chiral limit. The leading strong and 
electromagnetic contributions to the splittings dmx can 
be read off from the Lagrangians in Eq. (2.49). While 
both strong and electromagnetic terms contribute to the 
equidistant splitting 6m}, the non-equidistant one at this 
order of of pure electromagnetic origin. In Ref. (320), the 
values for ma and dmx” were determined from the most 
recent particle data group values for ma++ = 1230.80 + 
0.30 MeV and mao = 1233.45+0.35 MeV (324) together 
with the average mass ma = (ma++ + mMa+ + Mao + 
ma-)/4 = 1233 MeV from Ref. (325) which leads to 


dma = —5.342.0 MeV, dm\ = —1.7+2.7 MeV. 
(2.50) 
If the quark model relation (326) ma+—mMa,o = Mp—Mn 
is employed instead of using the average delta mass, the 
results change as follows 


dm_ = —3.9 MeV, dmA = —0.3+0.3 MeV. (2.51) 


which is consistent with Eq. (2.50). Notice that the val- 
ues for dmx” are of natural size. Indeed, based on naive 
dimensional analysis one expects |im,| ~ «M2/M, ~ 
8 MeV and |dm4| ~ e?M,/(47)? ~ 0.5 MeV. For a re- 
lated discussion on the delta mass splittings in chiral EFT 
with a particular emphasis on their quark mass depen- 
dence the reader is referred to Ref.(327). 

Having determined the values for the delta mass split- 
tings, it is a straightforward (but teddious) exercise to 
work out the leading A-contributions to the IB 2z- 
exchange potential. Notice that the A-contributions to 
the 17-exchange and contact potentials can be taken into 
account by a redefinition of various LECs and will, there- 
fore, not be discussed. The explicit expressions for the 
2n-exchange contributions can be found in Ref. (328). In 
Fig. 24 we show the CIB and CSB central, tensor and 
spin-spin potentials in coordinate space at N?LQ in the 
A-full theory in comparison with the N°L@ and N*L@ 
results in the A-less theory. While in the A-less theory, 
the leading and subleading class-II 27-exchange potential 
arises entirely from the pion mass difference M2, see the 
discussion in section II.B, in the A-full theory one also 
finds contributions proportional to 6m. Although these 
contributions are numerically small, they provide a clear 
manifestation of effects which go beyond the subleading 
order in the A-less theory. Furthermore, it is evident 
from Fig. 24 that the large portion of the N*L@ CIB 
27-exchange potential in the A-less theory is shifted to 
N°LQ in the theory with explicit A degrees of freedom 
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FIG. 24 Class-II (left panel) and class-III (right panel) 27- 
exchange potentials at N°L@ in the A-full theory (shaded 
bands) compared to the results in the A-less theory at 
N*°L@ (dashed lines) and N*L@ (dashed-dotted lines). The 
bands arise from the variation of dm and 6m according to 
Eq. (2.50). Notice further that the leading (ie. N°L@) contri- 
butions to Vis (r) and subleading (i.e. N*LQ) contributions 
to Vd" (r) vanish in the A-less theory. In all cases, the spectral 
function cutoff A = 700 MeV is used. 


leading to a more natural convergence pattern. Similarly, 
a comparison of the corresponding CSB (i.e. class-III) po- 
tentials in two theories also indicates towards a more nat- 
ural convergence in the A-full with the main part of the 
unnaturally large subleading contribution in the A-less 
theory being shifted to the lower order. Notice that the 
CSB 27-exchange potential at N3LQ in the A-full theory 
also receives contributions from the delta splitting 6m) 
which are still absent at N*+L@ in the A-less theory. For 
the central value, bm = —5.3 MeV, these contributions 
are numerically large and tend to cancel the ones driven 
by the nucleon mass difference leading to a significantly 
weaker resulting class-III 27-exchange potentials as com- 
pared to the ones at subleading order in the A-less theory. 
This can be viewed as an indication that certain higher- 
order IB contributions still missing at subleading order 
in the A-less theory are unnaturally large in the theory 
without explicit delta degrees of freedom. Last but not 
least, effects from virtual A-isobar excitation in the sub- 
leading my- and 27-exchange potential induced by addi- 
tional one-photon exchange are considered in Ref. (223). 

Inclusion of the A as an explicit degree of freedom has 
also important implications for the IB 3NF (320). As 
discussed in the previous section, the strongest IB 3NF 
arises from taking into account the charge-to-neutral pion 


mass difference in the 27-exchange 3N diagrams. In the 
A-less theory, the resulting charge-symmetry conserving 
3NF, being formally subleading (N*LQ), is enhanced by 
the large values of the LECs cz 4. In the A-full the- 
ory, the main part of this strong 27-exchange contribu- 
tion appears already at leading order (N°LQ) indicating 
a more natural convergence pattern. In addition to this 
obvious effect, one obtains further 27-exchange contribu- 
tions at N°LQ driven by the delta and nucleon mass split- 
tings 6m (charge-symmetry conserving) and dm), dm 
(charge-symmetry breaking). A close inspection of the 
resulting expressions, which are all proportional to A~?, 
reveals that they are exactly reproduced in the A-less 
theory by the saturation of the sub-subleading isospin- 
conserving (d;-terms in Lo) and IB zaNWN vertices. 
Numerically, the strengths of these CSB 3NFs due to in- 
termediate delta excitation turn out to be rather small, 
om — 36my|g2,h2,Me/(43202F4A2) ~ 3 keV, which, 
however, is comparable to the typical size of the remain- 
ing leading CSB 3NF (211), g4M4/(2567?F*) ~ 7 keV. 


E. Few-nucleon reactions involving pions 


Few-nucleon reactions involving pions such as 
eg. td — md (329), yd > 7w°d (330; 331), 7?He — 
mw He (332), ~d > ynn (333), yd — atnn (334) and 
NN — NNrz (335), where only some of the most recent 
references are given, provide another fascinating test- 
ing ground for the chiral EFT framework. The calcu- 
lations typically utilize the distorted-wave Born approx- 
imation using transition operators derived in chiral EFT 
and employing either phenomenological or chiral-EFT- 
based wave functions for the few-nucleon states following 
Weinberg’s original proposal (336). An important new 
ingredient in these applications is the appearance of the 
momentum scale p = \/mM,, associated with real pion 
production which has to be explicitly taken into account 
and requires an appropriate modification of the power 
counting (337-339). Such a modified ordering scheme 
was proposed in (283) and applied in (340) to calculate 
the pion production operator in NN collisions at thresh- 
old at NLO. Notice that the rather high energies and 
momenta of the nucleons in the initial state require the 
inclusion of the A-isobar as an explicit degree of free- 
dom. As a characteristic feature of the modified power 
counting scheme, one observes the appearance of half- 
integer powers of the small parameter y = M,/m in 
the expansion of the transition operators. One also finds 
that some pion loop contributions are promoted to signif- 
icantly lower orders compared to what is expected from 
Weinberg’s original power counting. An application of 
the modified power counting to P-wave pion production 
in NN collisions up to N?LO is carried out in Ref. (283). 
At this order, only tree diagrams have to be considered. 
The calculations showed a satisfactory agreement with 
the data and also demonstrate the feasibility to extract 
the LEC D which enters the leading 3NF, see Eq. (2.22), 
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FIG. 25 LO (dashed line) and NLO (solid line) results for the 
total cross section for the reaction pp — da* in comparison 
with the data from (342) (open circles), (343) (filled circles) 
and (344) (filled squares). The hatched area gives the esti- 
mated uncertainty at NLO. Figure courtesy of C. Hanhart. 


from this reaction. Notice, however, that concerns have 
been raised in Ref. (284) regarding the convergence of the 
chiral expansion in this reaction. For S-wave pion pro- 
duction, one-loop diagrams already start to contribute 
at NLO. As pointed out in Ref. (341), it is important to 
properly separate the truly irreducible contributions in 
the loop diagrams from the reducible ones in order for 
the resulting pion production operator to be renormaliz- 
able. Numerically, the NLO loop diagrams were found to 
provide an important contribution to the cross section for 
pp — dr+. Parametrizing the near-threshold cross sec- 
tion for this reaction as o = an+O(n*), where 7 denotes 
the outgoing pion momentum in units of its mass, the 
pion S-wave contribution at LO and NLO was found to 
be ab© = 131 ub and aN"© = 220 pb, respectively (341). 
The result at NLO agrees nicely with various existing 
data sets, see Fig. 25. A comprehensive review of meson 
production reactions in nucleon-nucleon collisions can be 
found in Ref. (345). 

Pion production reactions in few-nucleon systems also 
proved useful to study isospin-violating effects. Recent 
measurements of the forward-backward assymmetry in 
the process pn — dr (346) and the total cross section in 
the reaction dd + ar® (347) yielded a clear evidence of 
charge-symmetry breaking and serve as excellent testing 
ground to study isospin violation in the nuclear force and 
the corresponding transition operators. The first steps 
towards the theoretical understanding of these reactions 
have been taken in Refs. (348) and (349; 350), respec- 
tively. Notice that the appearance of the four-nucleon 
continuum states makes the theoretical analysis of the 
process dd — ar® particularly challenging. We further 
emphasize that new data on this reaction will be provided 
by WASA at COSY (351). Further details on these stud- 
ies and related issues can be found in a recent review 
article (352). 

The role of the momentum scale p = \/mM, and the 
related issue of the nucleon recoil effects in reactions such 


e.g. md scattering and pion photo- and electroproduction 
off the deuteron was investigated in the context of chiral 
EFT in Refs. (353-355). In particular, it was realized 
that the importance of the recoil effects in a given pro- 
cess is directly connected to the Pauli principle for the 
nucleons in the intermediate states. Notice further that 
the reaction yd — nna* (334) and the similar process 
md — ynn (333) were proposed as a tool to extract the 
value of the neutron-neutron S-wave scattering length. 
For more details on these and related reactions the reader 
is referred to Ref. (356). 


F. Hyperon-nucleon & hyperon-hyperon interactions 


The effects of strange quarks in nuclear matter can e.g. 
be tested through the determination of the properties of 
so-called hyper-nuclei, in which one (or two) nucleon(s) 
is (are) substituted by a hyperon (hyperons). Such 
hyper-nuclei are produced by strangeness-exchange re- 
actions, by associated strangeness production or by elec- 
troproduction reactions at many labs world-wide, such 
as CERN, BNL, KEK, DA®NE, JLab, MAMI and GSI, 
see e.g. Ref. (357). More generally, nuclear physics 
with strange quarks has a broad impact on contempo- 
rary physics since it lies at the intersection of nuclear 
and elementary particle physics. Moreover, it has signifi- 
cant implications to the astrophysics of compact objects. 
Recent progress in the field is reviewed in Ref. (358). 

The hyperon-nucleon (YN) interaction is at the heart 
of the hyper-nuclear binding and thus a precise deter- 
mination of its various components is of utmost impor- 
tance. Here, the situation is quite different compared 
to the two-nucleon case. The data base on YN scat- 
tering is quite poor, thus a partial wave analysis is not 
available and in any theoretical approach one must di- 
rectly compare to data. The poor status of our infor- 
mation on the YN interaction is most clearly reflected in 
the present knowledge of the AN scattering lengths. E.g 
Ref. (359) gives for the singlet (S) and triplet (T) scatter- 
ing lenghts a, = —1.8'7:5 fm, a, = —1.87)'3 fm, whereas 
in the six variants of the Nijmegen soft-core potential 
model a, varies between —2.5...—0.7fm and a; between 
—2.2...—1.8fm (360). In the most modern version of 
the Julich meson-exchange model one finds a, ~ —2.6 fm 
and a; ~ —1.7fm (361). However, for the EFT approach 
it is important to note that all these values are of natural 
size. For a proposal to extract these scattering lengths 
with high precision from production data, we refer to 
Ref. (362). Furthermore, since the masses of the A and 
the & hyperons are only about 75 MeV apart, the cou- 
pling between the AN and the UN channels needs to be 
taken into account. Moreover, for a sensible compari- 
son with experimental data, it is preferable to solve the 
scattering equation in the particle basis because then the 
Coulomb interaction in the charged channels can be in- 
corporated. 

The hyperon-nucleon YN interaction has not been in- 
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FIG. 26 One-pseudoscalar-meson-exchange diagrams at LO 
for the hyperon-nucleon interaction. 


vestigated using EFT as extensively as the NN interac- 
tion. Hyperon and nucleon mass shifts in nuclear mat- 
ter, using chiral perturbation theory, have been studied 
in (363). These authors used a chiral interaction contain- 
ing four-baryon contact terms and pseudoscalar-meson 
exchanges. The hypertriton (a bound state of a proton, 
a neutron and a A) and Ad scattering were investigated 
in the framework of an EFT with contact interactions 
(364). Korpa et al. (365) performed a next-to-leading or- 
der (NLO) EFT analysis of YN scattering and hyperon 
mass shifts in nuclear matter. Their tree-level ampli- 
tude contains four-baryon contact terms; pseudoscalar- 
meson exchanges were not considered explicitly, but 
SU(3) breaking by meson masses was modeled by incor- 
porating dimension two terms coming from one-pion ex- 
change. The full scattering amplitude was calculated us- 
ing the Kaplan-Savage-Wise resummation scheme. The 
hyperon-nucleon scattering data were described success- 
fully for laboratory momenta below 200 MeV, using 12 
free parameters. Some aspects of strong AN scattering 
in effective field theory and its relation to various formu- 
lations of lattice QCD are discussed in (50). 


Within the Weinberg counting scheme, a detailed in- 
vestigation of the YN interaction at LO was presented 
in (366). At LO, the YN potential is given by one- 
pseudoscalar-Goldstone-boson exchange diagrams, cf. 
Fig. 26, and contact interactons without derivatives. The 
spin-space part of the one-pseudoscalar-meson-exchange 
potential resulting from the LO SU(38) effective chiral 
meson-baryon Lagrangian is, 


V BiB. BiB, _ 


—fB, BP lB.ByP 


where k is the momentum transfer, P = 7, K,7, and the 


_Channel_Htsospinf er} _K by 


TABLE I The isospin factors for the various one— 
pseudoscalar-meson exchanges contributing to the hyperon- 
nucleon interaction. 


fp,B,P, fB.B,p are the appropriate coupling constants 


funx = f, fNNns = wl4a— Uf, 
fanxK = =a +20)f; fear = —U— 20)F, 
fEEns = = (1 +2a)f, fax = (4a = 1h, 
faonr = a Lays, fexns = All—alf, 
fenx = (1-2a)f, foun = 2af, 
faans = —Rll -a)f, fexx= —f. 

(2.53) 


in terms of the coupling constant f = g4/2F, and the 
F'/(F + D)-ratio a (367). The corresponding isospin fac- 
tors for the various channels multiplying the spin-space 
part of the potential in Eq. (2.52) are tabulated in Tab. I. 
It is important to stress that while the interaction po- 
tential at LO is SU(3) symmetric, the kinematics of the 
various channels and the masses of the exchanged mesons 
are to be taken at their physical values. The LO contact 
terms for the octet baryon-baryon interactions, that are 
Hermitian and invariant under Lorentz transformations, 
are given by the SU(3) invariants, 


LY = Ci (BB, (TB), TiB),) ; 
fe = GP BaV B) BiB) 5 
Le = C3 (B. (Ti B),) (By (TiB),) - (2.54) 


Here, a and b denote the Dirac indices of the particles, 
B is the usual irreducible octet representation of SU(3) 
given by 


x? A us 
Pa Lee ae P 
B= os ++, on ae ODS 
—E- m0 28 


and the brackets denote taking the trace in the three- 
dimensional flavor space. As an example, we display the 
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resulting partial wave potentials for AN > AN 


Vig = 4r 


| (C8 - 3Ch) + $ (C3 - 30%) 
2(c8 -3c3)| 
n|5 (Ch +04) + (C3 +C3) 


2(Cg+ c3)| ; (2.56) 
Similar expression for the isospin-1/2 and 3/2 UN — UN 
and the AN — XN potentials are given in Ref. (366). 
Note that only 5 of the {8} x {8} = {27} + {10} + 
{10*} + {8}, + {8}. + {1} representations are relevant 
for NN and YN interactions, since the {1} occurs only 
in the AA, =N and ©* channels. Equivalently, the 
six contact terms, C$, Cp, C2, C#, C3, C#, enter 
the NN and YN potentials in only 5 different combi- 
nations. These 5 contact terms need to be determined 
by a fit to the experimental data. The resulting chi- 
ral potential VY0 = Vopr + Veont in the Lippmann- 
Schwinger equation is regulated with a regulator function 
fa(p,p") = exp(—(p! + p“)/A4), where the cut-off A is 
varied between 550 and 700 MeV. A fit to 35 low-energy 
data (total cross sections from Refs.(359; 368-370) for 
Ap > Ap, 7p > Ap, +p — S*p and S7p > Yn 
with hyperon lab momenta between 110 and 300 MeV 
and the inelastic capture ratio at rest (371)) gives a good 
descripton of the data, see Fig. 27, with contact inter- 
actions of natural size. Note the strong cusp effect in 
Ap scattering at the opening of the ©*n threshold at 
Plab X 600 MeV (Fig. 27b). The chiral EFT also yields a 
correctly bound hypertriton (372) if one fixes the relative 
strength of the singlet and triplet S-waves accordingly. A 
Ap singlet scattering length of —1.9 fm leads to the cor- 
rect binding energy. The corresponding triplet scattering 
length is af? = —1.2 fm and in the +p system, one finds 
a, = —2.3 and a; = —0.7fm. It is quite astonishing that 
with six parameters only (5 LECs and the cut-off A) one 
achieves a quite satisfatory desciption of the admittedly 
not very preciese YN scattering data. Clearly, a NLO 
calculation should be performed and fits should be done 
simultaneously to YN and NN data. For a more detailed 
discussion of these results and a comparison to more con- 
ventional approaches to the YN interaction, we refer the 
reader to Refs. (366; 372). 

The experimental situation on baryon-baryon scatter- 
ing with S = —2, ie in the YY and the EN chan- 
nels, is even poorer. Only very recently doubly strange 
baryon-baryon scattering data at lower energies, below 
Plab = 0.8 GeV, were deduced for the first time (373; 374). 
An upper limit of 24 mb at 90% confidence level was pro- 
vided for elastic =~ p scattering, and for the =" p — AA 
cross section at Piab = 500 MeV a value of 4.3*§3 mb 
was reported (374). Within LO chiral EFT, baryon- 
baryon scattering was analyzed by Polinder et al. (375). 
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FIG. 27 Total cross sections as a function of piab for Ap — Ap 
and *p — “~p. The shaded band is the LO EFT result for 
A = 550...700 MeV, the dashed curve is the Jiilich 04 model 
(361), and the solid curve is the Nijmegen NSC97f model 
(360). Note that the total cross sections for U*p — U*p 
are obtained by integrating the differential data in a limited 
angular range, see Ref.(369). 


The contact terms and the couplings of the pseudoscalar 
mesons to the baryons are related via SU(3) symmetry 
to the S$ = —1 hyperon-nucleon channels. There is one 
additional contact interaction whose strength was var- 
ied within natural bounds in the AA — AA channel. 
This fixes its contribution in all other S = —2 baryon- 
baryon channels because of SU(3) symmetry. As a con- 
straint, the information deduced from the recent candi- 
date for ,’He with a low binding energy (376), the so- 
called Nagara event, that suggests that the AA interac- 
tion should be only moderately attractive, was imposed. 
For a fixed cut-off, the prediction of the =~p — = 7p 
and the =-p — AA total cross section in comparison 
to the available data are shown in Fig. 28. The re- 
sulting AA scattering length in the 1S channel is a, = 
—1.83...— 1.38fm. For comparison, in the Nijmegen 
ESC04 model one finds in this channel a, = —1.32fm 
(377) and in the constituent quark model of Fujiwara et 
al. (378) one has a, = —0.81fm. Note also that this 
contact interaction does not contribute to certain chan- 
nels, so that at LO one can make parameter-free predic- 
tions for N+Ht 3 UY, By > 2% and B°p > UTA. 
It is expected that in the coming years better-quality 
data on the fundamental =N and YY interactions as well 
as much more information about the physics of hyper- 


33 


o (mb) 
8 
T 
L 


0 1 1 L 
0 200 400 600 800 


Plab (MeV/c) 


o (mb) 


0 1 1 i 
0 200 400 600 800 


Plab (MeV/c) 


FIG. 28 Total cross sections as a function of pap for =~ p > 
=~ p (upper panel) and =p — AA. The shaded band is the 
LO EFT result for A = 600 MeV and varying the LEC Ci) 
of the additional singlet contact term within natural bounds 
gievn a mildly attractive interaction in the ' So channel of the 
AA interaction. 


nuclei will become available at the new facilities J-PARC 
(Japan) and FAIR (Germany). The chiral EFT devel- 
oped in Ref. (375) can then be used to analyze these 
upcoming data in a model-independent way. 


G. Nuclear lattice simulations 


Once the chiral nuclear forces are determined and the 
low-energy constants appearing in the nuclear forces are 
fitted (in the two and three-nucleon sector) one can make 
predictions in the four- and more-nucleon sectors based 
on chiral EFT. However, the explicit numerical treatment 
of e.g. the Yakubowsky equations for more than four 
nucleons is a very difficult task. One possible scheme 
to solve the many-body problem is to put the chiral ef- 
fective potential on the lattice and apply the powerful 
Monte-Carlo techniques which are already developed to 
high degree in lattice QCD. One unique feature of the lat- 
tice effective field theory approach is the ability to study 
in one formalism both few- and many-body systems as 
well as zero- and nonzero-temperature phenomena. A 
large portion of the nuclear phase diagram can be stud- 
ied using exactly the same lattice action with exactly the 
same operator coefficients. A second feature is the com- 
putational advantage of many efficient Euclidean lattice 
methods developed for lattice QCD and condensed mat- 
ter applications. This includes the use of Markov chain 


Monte Carlo techniques, Hubbard-Stratonovich transfor- 
mations, and non-local updating schemes such as a hy- 
brid Monte Carlo. A third feature is the close theoretical 
link between nuclear lattice simulations and chiral effec- 
tive field theory. One can write down the lattice Feyn- 
man rules and calculate lattice Feynman diagrams using 
precisely the same action used in the non-perturbative 
simulation. Since the lattice formalism is based on chi- 
ral effective field theory we have a systematic power 
counting expansion, an a priori estimate of errors for 
low-energy scattering, and a clear theoretical connection 
to the underlying symmetries of QCD. The first studies 
combining lattice methods with effective field theory for 
low-energy nuclear physics looked at infinite nuclear and 
neutron matter at nonzero density and temperature (see 
Sec. III.E). Most of the formalism for chiral EFT on the 
lattice was developed in (379). Nuclear lattice simula- 
tions were used to study the triton at leading-order in 
pionless effective field theory with three-nucleon interac- 
tions (380). 

Let us briefly discuss the principles underlying such 
nuclear lattice simulations (for a detailed discussion, we 
refer to Refs. (379; 381)). In this framework, nucleons 
are represented as point-like Grassman-fields and pions 
as point-like instantaneous pseudoscalar fields. The lat- 
tice is defined by a volume L? x L;, with L(L,) the 
spatial (temporal) size. The corresponding lattice spac- 
ings are called a and az, respectively. Typically, calcu- 
lations are carried out using a lattice length DL ~ 20 fm 
and the lattice spacing a ~ 2 fm which corresponds to 
the cutoff A = a1/a ~ 300MeV. At present, computa- 
tional resources prevent one from using smaller lattice 
sizes (larger UV cutoffs). Still, the use of various forms 
of improved actions allows one to access the systematic 
errors inherent in such simulations. 

The basic quantity in nuclear lattice simulations is the 
correlation function. For A nucleons in Euclidean space 
it is defined by 


Za(t) — 


where |W 4) refers to a Slater determinants for A free nu- 
cleons, H is the Hamiltonian of the system and 7 the Eu- 
clidean time. The ground state energy of the A-nucleon 
system can be derived from the asymptotic behavior of 
the correlation function for large 7, 


(WU 4| exp(—TH)|W a), (2.57) 
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The expectation value of any normal ordered operator O 
can be derived in a similar way by 


23 (7) 


(W2)O|WA) = him Zar)” (2.59) 
ZQ(t) = (Wal exp(—7TH/2)O exp(—TH/2)|WVa), 


where |W%) denote the ground states of A-nucleon sys- 
tem. It is convenient to describe NN contact interac- 
tions by standard bilinear nucleon density operators us- 
ing the Hubbard-Stratonovich transformation. Using the 
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relation exp(p?/2) ~ f ds exp(—s?/2 — sp) one can ex- 
press terms quadratic in the nucleon density operator p as 
terms linear in p in the presence of auxiliary background 
fields (collectively denoted by s). In this representation, 
the full correlation function is related to the path integral 
over pions and auxiliary fields, 


ZA ) =] Ds II Dr; Ds; exp(—Saa — Sss) 


I=1,2,3 


x (Wy4|M“t-) (ay, 8, 87) --- MO (az, 8, 87)|Va) « 


(2.60) 


Here S,, and S, are free actions for pions and aux- 
iliary fields s,s; (where s/s; couples to the isospin- 
independent/dependent nucleon bilinear), J denotes 
isospin indices and WV is an (irrelevant) normalization 
constant. M ”) is a transfer matrix defined as an n’th 
step in the temporal direction. We remark that the am- 
plitude (W4|M*-) (ry, 8;) +» MO (az, 5;)|Wa) is just 
a Slater determinant of single nucleon matrix elements 
Mi; with 4,9 = 1,...,A. 

To be specific, we give here the leading order action 
starting with the free theory. The presentation here is 
somewhat sketchy. For an extensive discussion see (381). 
The free actions for the auxiliary fields and the pions are 


3 
a Daa 
yom 


Sss(5, 81) _ oe (71) 


San (tr) = (—A + M2)r7(a Oe 


(2.61) 


where M, is the physical pion mass and a, = a;/a. 
For nucleons one may use an O(a*) improved free lat- 
tice Hamiltonian defined by 


ive 4 Shea 
eee = oe S- fr [al (fs) (aig (As + kls) 


k=0 #05 ,4,7 


+4;,; (iis — kf,))], 
where the operators al, j(vis) and aij (7is) are the nucleon 
creation and annikilation operators, 7, are spatial coor- 
dinates, i, are spatial unit vectors, the indices 7 and 7 rep- 
resent spin and isospin indices, respectively, and the co- 
efficients f, are fo,1,2,3 = 49/2, —3/4, 3/40, -1/180. To 
define the interactions one introduces nucleon-density op- 
erators with different spin/isospin polarizations 


p* “(hs) = 2 ahlfaaes (fe) 


(2.62) 


po ia oe a; Jez] sr.405.5Gi5)5 (2.63) 
45,5" 

al aso = =" 

Pr.§ ( )= S- al, 4, (fis)[os]v alr] j4i,5 (Ts). 


Pspeue! Sp 
U,V 5I5d 


The transfer matrix for n;-th step has, besides the free 
part, two important contributions: 


M(™) _ exp ~ Hist 
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Here :: denotes normal ordering. The first long-range 
contribution ~ gy includes the instantaneous pion- 
nucleon interaction and describes the one-pion-exchange 
in the leading-order effective potential. The second short- 
range contribution corresponds to the NN contact in- 
teractions. The low-energy constants C = C's — 2C'r 
and C; = —Cr (cf. sec. IIB) have different signs, 
C <0,C; > 0. With these signs the pion-less theory 
can be shown to have no sign-oscillations if the number 
of protons and neutrons are equal and they stay pair-wise 
in isospin-singlet states. In this case the multiplication 
with 72 of the single-nucleon matrix elements M from left 
and right is well defined and gives 72Mt2 = M*. For this 
reason, the determinant of M is real, det M* = det M. 
Since Ty is antisymmetric, the eigenvalues of M are dou- 
bly degenerate. This leads to a positive Slater determi- 
nant (382; 383) 


det M > 0. (2.65) 


The introduction of pions causes small sign oscillations 
which, however, are not severe and appear to be sup- 
pressed. 

To perform simulations in a most efficient way one ex- 
ploits the approximate SU(4)-Wigner (384) symmetry in 
the NN system. The symmetry transformation is given 
by independent rotations of the spin and isospin degrees 
of freedom. 


ON =ayo't’N with oM =(1,¢), 7% =(1,7). 

(2.66) 
One can show that in the limit where the NN S-wave scat- 
tering lengths approach infinity the two-nucleon system 
becomes invariant under the SU(4)-transformation (385). 
The SU(4)-breaking corrections come from the finite 
the scattering length and higher order terms in the 
chiral expansion, these are of order O((1/a(?Si) — 
1/a(*So)),q/Ay). Since the NN scattering lengths are 
very large, the SU(4)-breaking corrections appear to be 
small. This fact can be used to improve the performance 
of the lattice simulations. The SU(4) symmetric transfer 
matrix is given by 


M) = :exp | Hirxc 


+ V/—Caz S~ slits, nz)p"“(iis)] + . (2.67) 
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FIG. 29 Overview of the various pieces of the transfer matrix 
calculation. 


In this case there are no sign oscillations for an even 
number of nucleons (386) and one only has one auxiliary 
field such that the simulations become much cheaper. 
Although there is no positivity theorem for odd numbers 
of nucleons, sign oscillations seem to be suppressed also 
in systems with odd number of nucleons because it is 
only one particle away from an even system with no sign 
oscillations. Since the final result is close to the one pro- 
duced by a SU(4)-symmetric simulation it pays to divide 
a simulations into three parts. To simulate the expecta- 
tion value of some observable one uses SU(4)-symmetric 
transfer matrices in the first and the last Ly, steps in 
order to filter the low-energy signal. After this filtering, 
one starts the simulation with the complete (realistic) 
transfer matrices. A schematic overview of the transfer 
matrix calculation is shown in Fig. 29. Having set up 
the transfer matrix, one utilizes the hybrid Monte-Carlo 
(HMC) method (387) to update the field configurations. 
More specifically, one introduces the conjugate fields p,.,, 
Ds, Ps, and uses molecular dynamics trajectories to gen- 
erate new configurations for the fields pz,, Pps, Ps;, TI, 3; 
sr which keep the HMC Hamiltonian 


Hnsc = 5 (Db. (a) + 93,0] +230)) 
nm I 


+ V(a1, 8,81), (2.68) 
constant, where the HMC potential is defined by 
V (m1, 8, $1) = Saar + Sss — log| det M|. (2.69) 


Upon completion of each molecular dynamics trajec- 
tory, an Metropolis accept or reject step for the new 
configuration according to the probability distribution 
exp(—Humc) is applied. This process of molecular dy- 
namics trajectory and Metropolis step is repeated many 
times. 

Already at LO promising results for binding energies, 
radii and density correlations for the deuteron, triton and 
“He are obtained (381). On a 5° lattice, the triton bind- 
ing energy agrees with experiment within 5% and the 
triton root-mean-square radius is accurate to 30%. The 
binding energy for *He is within 25% of the experimen- 
tal value while the root-mean-square radius agrees within 
10%. Note, however, that one has to overcome a zero- 
range clustering instability that appears for 4 (or more 
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FIG. 30 NN S-wave phase shifts and mixing angles versus 
center-of-mass momentum with actions LO; and LO2. 


particles) which is mostly a combinatorial effect when 
more than two particles occupy the same lattice site (for 
studies of this in other systems, see (388) and references 
therein). To overcome this problem, one can e.g. smear 
out the contact interactions with a Gaussian. Such terms 
improve the lattice action and are formally of higher or- 
der (for more details, see e.g. Ref. (381). At LO, one can 
also study the feasibility of simulations for light nuclei 
with more than four nucleons. It was observed that for 
A < 10 the CPU time scales approximately linearly with 
the nucleon number A. 

At NLO there appear 9 LECs which can be fitted to 
the Nijmegen NN scattering data and deuteron proper- 
ties. Elastic scattering phase-shifts on the lattice are re- 
lated by Liischer’s formula to the energy levels of two- 
body states in a finite large volume cubic box with pe- 
riodic boundary conditions (see Sec. I.E). While this 
method is very useful at low momenta, it is not so use- 
ful for determining phase shifts on the lattice at higher 
energies and higher orbital angular momenta. Further- 
more, spin-orbit coupling and partial-wave mixing are 
difficult to measure accurately using Liischer’s method 
due to multiple-scattering artifacts produced by the pe- 
riodic cubic boundary conditions. In Ref. (389) a more 
robust approach to measure phase shifts for two nonrela- 
tivistic point particles on the lattice using a spherical wall 
boundary was proposed. The basic idea is to impose a 
hard spherical wall boundary on the relative separation 
between the two interacting particles at some chosen ra- 
dius. The reason for this spherical wall is to remove 
copies of the two-particle interactions due to the peri- 
odic boundaries on the lattice. This additional bound- 
ary condition allows for a direct extraction of the phase- 
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shifts and mixing angles from the finite-volume spectrum. 
For more details, see (389). Using the spherical wall 
method the values of 9 LECs were determined by match- 
ing three S-wave and four P-wave scattering data points, 
as well as deuteron binding energy and quadrupole mo- 
ment. In Fig. 30 the NN S-wave phase-shifts and the 
39, —3 Di mixing angle e€ for two different actions, called 
LO, and LOg, are displayed. The action LO, is the one 
given in Eq. (2.64). In the action LOz2 the contact inter- 
actions are smeared by a Gaussian. The two actions are 
identical at leading order and differ only by higher-order 
terms, thus given an estimate of the higher order correc- 
tions. As can be seen from Fig. 30, the results of the lat- 
tice simulations are in a good agreement with the partial 
wave results for momenta smaller than ~ 100 MeV. Devi- 
ations between the two results for different actions appear 
merely at larger momenta and are consistent with the 
expected higher order effects. In the mean time, a novel 
action with spin-isopsin projected smearing has been de- 
veloped that gives a good description of the partial waves 
up to momenta of the order of the pion mass. 

At N?LO three-body forces start to show up which de- 
pend on two constants. These LECs can be determined 
from a fit to neutron-deuteron scattering data in the spin- 
1/2 doublet channel and the triton binding energy. These 
simulations show a very natural convergence pattern with 
increasing chiral order. For a box length ~ 15 fm the 
volume dependence already becomes very small and the 
binding energy approaches its physical value. This is 
consistent with our expectation that the volume depen- 
dence in nuclear lattice EFT simulations should become 
weak for L ~ 20 fm. In Fig. 31 the S-wave phase-shifts in 
the spin-3/2 quartet channel versus the square of relative 
momentum are shown. This channel was not taken into 
account in the fit procedure. Again one observes a very 
nice convergence with increasing chiral order. The pre- 
dictions are located between the experimental data for 
proton-deuteron and neutron-deuteron scattering data. 
Since isospin-breaking was not taken into account in the 
simulations the results are very satisfactory. At the same 
order, the (Coulomb-corrected) binding energy for *He is 
overpredicted by 5%, which is consistent with the ex- 
pected theoretical accuracy of these simulations. 

The results of these studies demonstrate that lattice 
EFT is a promising tool for a quantitative description 
of light nuclei. In the future, it is planned to perform 
N?LO Monte-Carlo simulations of light nuclei and probe 
neutron matter with larger number of neutrons in a box, 
see also Sec. HI.E. In addition, more detailed studies of 
finite size effects and further improvements of the lattice 
action are called for. 


H. Quark mass dependence of nuclear forces and IR limit 
cycle in QCD 


The quark mass dependence of the chiral NN inter- 
action was calculated to next-to-leading order (NLO) 
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FIG. 31 Neutron-deuteron scattering S-wave phase-shifts in 
the spin-3/2 quartet channel versus the square of relative mo- 
mentum. The data for proton-deuteron and neutron-deuteron 
scattering are taken from Ref. (390). 


in the chiral counting in Refs. (154; 391). At this or- 
der, the quark mass dependence is synonymous to the 
pion mass dependence because of the Gell-Mann-Oakes- 
Renner relation: M? = —(my + ma)(0|au|0)/F? , where 
(O|tiu|0) + (—290 MeV)? is the quark condensate. In the 
following, we will therefore refer only to the pion mass 
dependence which is more convenient for nuclear applica- 
tions and treat the pion mass as a parameter that can be 
varied by adjusting the values of the quark masses. In the 
work of Refs. (154; 266; 391), the pion mass dependence 
of the nucleon-nucleon scattering lengths in the 3S;-3D, 
and ‘So channels as well as the deuteron binding energy 
are calculated. To next-to-leading order (NLO) in the 
chiral power counting, the NN potential can be written 
as 

VxLo a YOPE Ay yVTPE ze Voont , (2.70) 
where VOPE, YTPE. and Vt refer to the one-pion 
exchange, two-pion exchange, and contact potentials, 
respectively. Explicit expressions can be found in 
Ref. (154) and section II.B. 

In principle, the pion mass dependence of the chiral 
NN potential is determined uniquely. However, the ex- 
trapolation away from the physical pion mass generates 
errors. The dominating source are the constants Csr 
and Ds r in V°°" which give the corrections to the LO 
contact terms « M? and cannot be determined indepen- 
dently from fits to data at the physical pion mass. A 
smaller effect is due to the error in the LEC dj6, which 
governs the pion mass dependence of g4, from the chiral 
pion-nucleon Lagrangian which is enhanced as one moves 
away from the physical pion mass. Both effects generate 
increasing uncertainties as one extrapolates away from 
the physical point. 

In the calculation of Ref. (154), the central value 
dig = —1.23 GeV~? was used which is the average of 
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FIG. 32 Deuteron binding energy as a function of the pion 
mass M,. The shaded areas correspond to the allowed val- 
ues. The light shaded band gives the uncertainty due to the 
unknown value of the LECs Dg,r using the central value 
dig — 1.23 GeV~?. The dark shaded band gives the uncer- 
tainty if, in addition to the variation of Ds,r, the LEC die 
is varied in the range from dig — 0.91 GeV~? to dig — 1.76 
GeV~? given in Ref. (392). The heavy dot shows the binding 
energy for the physical value of the pion mass. 


three values given in Ref. (392). In addition, a variation 
of the LEC dg in the range dig = —0.91...—1.76 GeV~? 
as given in Ref. (392) was employed. The size of the two 
constants Dg and Dr was constrained from naturalness 
arguments. It was argued that the corresponding dimen- 
sionless constants F2A> Ds,r can be expected to satisfy 
the bounds: 

-3< F2M Der <3, (2.71) 
where A, ~ 1 GeV is the chiral symmetry breaking 
scale. A more conservative error estimation was given 
in Ref. (393). We note that Refs. (266; 391) allowed 
for a larger variation of these LECs. However, the 
bounds (2.71) are in agreement with resonance satura- 
tion estimates and similar relations are obeyed by the 
LECs whose values are known (280). For the con- 
stants Csr, e.g., this leads to Cs = —120.8 GeV~? and 
Cr = 1.8 GeV~? corresponding to the dimensionless co- 
efficients F?Cs = —1.03 and F2Cr = 0.02, respectively. 
The unnaturally small value of F2C’p is a consequence of 
the approximate Wigner SU(4) symmetry. 

The ranges from Eq. (2.71) were used to estimate the 
extrapolation errors of two-nucleon observables like the 
deuteron binding energy and the spin-singlet and spin- 
triplet scattering lengths in Ref. (154). The resulting 
pion mass dependence of the deuteron binding energy is 
shown in Fig. 32. In the chiral limit the binding energy is 
of natural size, Bp ~ F?/m ~ 10 MeV. Note, however, 
that in the calculation of (266; 391) the assumed larger 
uncertainties in the LECs prevent one from making a def- 
inite statement about the binding of the deuteron in the 
chiral limit. For pion masses above the physical value 
the differences between the two calculations are consid- 
erably smaller. The recent study of Mondejar and Soto 
seems to indicate that two-loop diagrams generate a pe- 
culiar quark mass dependence of the contact interactions 
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FIG. 33 Inverse of the S-wave scattering lengths in the spin- 
triplet and spin-singlet nucleon-nucleon channels as a function 
of the pion mass M,. Filled triangles and rectangles show the 
lattice calculations from Refs. (42; 43) and (49), respectively. 
For remaining notation see Fig. 32. 


which are parametrically large (138). The influence of 
these effects on the quark mass dependence of, e.g., the 
deuteron binding energy remains to be worked out in de- 
tail. In Fig. 33, we show the inverse scattering lengths in 
the spin-triplet and spin-singlet channels from Ref. (154) 
together with some recent lattice results (49). However, 
the errors and pion masses are still too large to draw any 
conclusions about the physical point. 

Fig. 33 also shows that a scenario where both inverse 
scattering lengths vanish simultaneously at a critical pion 
mass of about 200 MeV is possible. For pion masses be- 
low the critical value, the spin-triplet scattering length 
would be positive and the deuteron would be bound. As 
the inverse spin-triplet scattering length decreases, the 
deuteron would becomes more and more shallow and fi- 
nally would become unbound at the critical mass. Above 
the critical pion mass the deuteron would exist as a shal- 
low virtual state. In the spin-singlet channel, the situ- 
ation is reversed: the “spin-singlet deuteron” would be 
a virtual state below the critical pion mass and would 
become bound above. It is unlikely that this scenario of 
both inverse scattering lengths vanishing simultaneously 
is realized in QCD at the physical values of the quark 
masses. However, based on this behavior it was conjec- 
tured that one should be able to reach the critical point 
by varying the up- and down-quark masses m, and mq 
independently because the spin-triplet and spin-singlet 
channels have different isospin (82). In this case, the tri- 
ton would display the Efimov effect which corresponds to 
the occurence of an infrared limit cycle in QCD. It is ev- 
ident that a complete investigation of this issue requires 
the inclusion of isospin breaking corrections and there- 
fore higher orders in the chiral EFT. However, a number 
of studies have investigated the universal properties of 
the limit cycle by considering specific values of Dg and 
Dr. 

In the exploratory study of the three-nucleon system 
(82), the mean values of the error bands from (154) were 
used as input for the three-body calculations in the pi- 
onless EFT. Even though both scattering lengths were 
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FIG. 34 Binding energies B3 of the triton ground and first 
two excited states as a function of M,. The circles, squares, 
and diamonds give the chiral EFT result, while the solid lines 
are calculations in the pionless theory. The vertical dotted 
line indicates the critical pion mass M¢"™ and the dashed 
lines are the bound state thresholds. 


large for the mean values, they did not become infinite 
at the same value of the pion mass and there was no ex- 
act limit cycle for this choice of parameters. However, 
different sets of values for Dg and Dr that lie within the 
bound given by Eq. (2.71) and cause the spin-singlet and 
spin-triplet scattering lengths to become infinite at the 
same value of the pion mass can be found. 


In Ref. (394), the properties of the triton around the 
critical pion mass were studied for one particular solu- 
tion with a critical pion mass M°"’ = 197.8577 MeV. 
From the solution of the Faddeev equations, the bind- 
ing energies of the triton and the first two excited states 
in the vicinity of the limit cycle were calculated for this 
scenario in chiral EFT. The binding energies are given in 
Fig. 34 by the circles (ground state), squares (first excited 
state), and diamonds (second excited state). The dashed 
lines indicate the neutron-deuteron (M, < Ms") and 
neutron-spin-singlet-deuteron (M, > M¢") thresholds 
where the three-body states become unstable. Directly 
at the critical mass, these thresholds coincide with the 
three-body threshold and the triton has infinitely many 
excited states. The solid lines are leading order calcula- 
tions in the pionless theory using the pion mass depen- 
dence of the nucleon-nucleon scattering lengths and one 
triton state from chiral EFT as input. The chiral EFT 
results for the other triton states in the critical region are 
reproduced very well. The binding energy of the triton 
ground state varies only weakly over the whole range of 
pion masses and is about one half of the physical value 
at the critical point. The excited states are influenced by 
the thresholds and vary much more strongly. 


These studies were extended to N?LO in the pion- 
less EFT and neutron-deuteron scattering observables in 
Ref. (395). It was demonstrated that the higher order 
corrections in the vicinity of the critical pion mass are 
small. This is illustrated in Fig. 35, where we show the 
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FIG. 35 Doublet neutron-deuteron scattering length res in 
the critical region computed in the pionless EFT. The solid 
line gives the LO result, while the crosses and circles show the 
NLO and N?LO results. The dotted lines indicate the pion 
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doublet scattering length ae in the critical region. The 


solid line gives the LO result, while the crosses and cir- 
cles show the NLO and N?LO results. The dotted lines 
indicate the pion masses at which a. i diverges because 
the second and third excited states of the triton appear 
at the neutron-deuteron threshold. These singularities in 
aoe (M,,) are a clear signature that the limit cycle is ap- 
proached in the critical region. Considering the specific 
case of the doublet scattering length at MM, = 190 MeV, 
we have 


a.!? (My = 190 MeV) = (93.18+-0.80+0.14) fm . (2.72) 


At this pion mass we are fairly far away from any singular 
points in the function at!? (Mz). Consequently, the pion- 
less results follow a natural convergence pattern with the 
expansion parameter yro, which is + 0.08 at this value 
of M,. For three-body scattering observables in the crit- 
ical region no calculation in chiral EFT is available. The 
above example shows that both approaches complement 
each other. If the pion mass dependence of one three- 
body observable is known, the pionless theory can be 
used to predict all other observables with high precision 
and less computational effort. Figure 35 demonstrates 
clearly that the pionless theory converges rapidly in the 
critical region. 


A final answer on the question of whether an infrared 
limit cycle can be realized in QCD can only be given by 
solving QCD directly. In particular, it would be very 
interesting to know whether this can be achieved by ap- 
proriately tuning the quark masses in a Lattice QCD 
simulation (396). Finally, we note that the dependence 
of nuclear binding on hadronic mass variation has also 
been studied based on the Argonne potential (397). 
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Ill. TOWARDS A MANY-BODY EFT FOR NUCLEI 


The EFT approach to the nuclear many-body problem 
is much less established than the one for the forces and 
few-nucleon systems. This is, on one hand, related to the 
apperance of new scales like the Fermi momentum or in- 
duced by collective excitations, and on the other hand to 
the computational problems related to solve the many- 
body problem. Thus, a variety of pathways are being ex- 
plored and here we can only give a brief overview about 
the existing attempts and their status. For a pedagog- 
ical review on the application of EFTs to finite density 
systems, we refer to Furnstahl et al. (398). 


A. In-medium chiral perturbation theory 


Early attempts to formulate in-medium chiral per- 
turbation theory where mostly triggered by the pio- 
neering paper of Kaplan and Nelson on kaon condensa- 
tion (399). Most of these calculations were based on chi- 
ral Lagrangians at most bilinear in the nucleon fields and 
performed in the mean-field approximation, NDN — 
PplrDi1 + pnTrD22 with pp(pn) the proton (neutron) 
density, D represents a generic differential operator in- 
cluding the coupling to pions and external sources, the 
trace runs over spinor indices and the subscripts run in 
flavor space. Proceeding in this way, one keeps track 
about the vacuum CHPT Lagrangians, but the chiral 
counting in the medium is lost, as e.g. nucleon corre- 
lations are not considered. The most elegant formulation 
of this approach based on the path-integral formulation 
is due to Wirzba and collaboraotrs, see e.g. (400-402). 

To go beyond the mean-field approximation, the in- 
medium genrating functional for pions coupled to nucle- 
ons and external sources was developed in Refs. (403; 
404). Leaving out multi-nucleon interactions, as system- 
atic in-medium CHPT can be developed by expanding 
around the nuclear matter ground state at asymptotic 
times and intergrating out the nucleon fields in the path 
integral representation, giving rise to the in-medium gen- 
erating functional (for the detailed derivation, we refer to 
Ref. (403)) 
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FIG. 36 Generalized in-medium vertices of lowest order. The 
thick solid lines correspond to insertions of a Fermi-sea and 
each circle to the insertion of the operator I as defined in the 
text. 
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where the operator A is defined as the difference between 
the full and the free Dirac operators, 


Lon = Ha) DCa)¥le) = Ha) Dole) — ACV)» 
with Do = iy,0"% — m, while the diagonal flavor ma- 
trix n(p) = diag(0(k\?) — |p|, 0(&k&? — |p|) parameterizes 
the upper cut-off of the three-momentum integrations in 
terms of the proton and neutron Fermi momenta, re- 
spectively. Furthermore J, is the unit operator in four 
dimension, E(p) the on-shell energy of a nucleon with 
mass m and v,a,8,p are vector, axial-vector, scalar and 
pseudoscalar sources. The resulting in-medium effective 
Lagrangian £ is given in terms of pions and external 
sources only and thus the problem is reduced to that of 
vacuum CHPT, with the important difference that the £ 
is non-covariant as well as non-local (for a general anal- 
ysis of the structure of non-relativistic but local EFTs 
see (405)). In particular, we note the appearance of the 
non-local vacuum vertex T = —iA(I, — Dj'A)~! that 
generates a geometric series in terms of the local inter- 
action operator A and the free Dirac propagator, with 
A itself being subject to the standard chiral expansion, 
A = AM + AQ) +..., see (145). The generalized in- 
medium vertices, cf Fig. 36, consist of several non-local 
vacuum vertices [ connected through the exchange of 
on-shell Fermi-sea states. These are the building blocks 
for the systematic expansion in small momenta, count- 
ing the Fermi momentum kr ~ 2M, at nuclear satura- 
tion as O(p). The in-medium chiral counting including 
the contributions from nucleon propagators can now be 
given. The choice of the counting scheme depends on 
the energy flowing through the nucleon lines, inducing a 
separate consideration of the so-called standard amd non- 
stantard cases. Let us consider first the former. Here, the 
energy flow is of order M, ~ O(p) and thus the nucleon 
propagator counts as Dj‘ ~ O(p~'). The chiral dimen- 
sion for a many-particle diagram with L, pion loops and 
Vr vacuum and/or in-medium vertices of dimension 6; is 


Vr 
vy =2L,+2+ 5 (5; —2). 


j=1 


(3.3) 
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Consequently, the lowest order in-medium contributions 
arise at O(p*) since the lowest order in-medium vertices 
have dimesnion four due to the four-momentum Dirac 
delta-function attached to any T vertex. The first correc- 
tions at NLO arise at O(p°), which should be contraested 
to the vacuum case where LO (NLO) is O(p?) (O(p*)). 
In the absence of multi-nucleon interactions, the break- 
down scale is A = V/67F,, ~ 700MeV for S-waves and 
A= V6rF, /ga =~ 560 MeV for P-waves. However, there 
is one subtlety with this power counting. Quite similar 
to what happens in case of the TPE NN interaction, the 
energy flowing into a nucleon line can vanish, so that the 
nucleon propagator scales as O(p~”). To deal with this 
non-standard case, one has to separately count the num- 
ber of nucleon lines with energy E < k#/2m and the 
normal lines with F ~ M,. The explicit expression for 
the modified counting index v can be found in (404). In 
this case, the breakdown scale is 617F?/2m ~ 270 MeV 
for S-waves and further reduced by a factor of 1/g% for 
P-wave interactions. Note that the so-defined in-medium 
CHPT not only encompasses but also transcends the so- 
called low-energy theorems of Refs. (406-408). 

We now discuss some results obtained in this scheme. 
The density dependence of the light quark condendates 
is given at NLO by 


_ _ 20. 4e5 
(Olu) = (te [1 = serra + Sp] 

2 = 20. Acs 
(Qdad]O) = (da) vec f\- 0 - SA] (3.4) 


with |Q) the nuclear matter background, p = (pp+pn)/2, 
P = (Pp — Pn) /2 are the isospin symmetric and asymmet- 
ric combinations of the proton and the neutron densities, 
while o is the pion-nucleon sigma-term, ¢ = —4c, M2? at 
O(p). The small isospin-breaking contribution is given 
in terms of the LEC cs = —(0.09 + 0.01) GeV. Further- 
more, the subscript vac refers to the vacuum value of the 
corresponding quantity. Higher order corrections in the 
density will be discussed in the next subsection. 

The propagation of pions in the medium can be an- 
alyzed by calculating the spectral relations between the 
energy w and the three-momentum q for on-shell neu- 
tral and charged pions. For symmetric nuclear matter 
(with density #) and in the chiral limit one obtains the 
dispersion law 


4b 
2 Pp 

Ww? = (1 = Fe ca) (3.5) 
Since the in-medium pion velocity 6 = dw/d|¢| = 


1 — 2fco/F? must be smaller than the velcoty of light 
(405; 409), this imposes the constraint cp > 0, which is 
satisfied by the actual value of this LEC. It was also es- 
tablished in (404) that for standard values of the LEC 
c3, chiral symmetry can account for the observed mass 
shift of the negatively charged pion in deeply bound pi- 
onic states in 7°’Pb (410; 411). At NLO, one obtains 


AM, =18+5MeV which is compatible with the exper- 
imental result, AM, = 23...27MeV. Of interest is also 
splitting of the temporal and space-like components of 
the pion decay constant, which read at NLO in symmet- 
ric nuclear matter 


F,=F, {1 ~£(0.26+ 0.04} 
Po 


F, = Fy {1 = f (1.93 + 007} (3.6) 


Po 

with po the nuclear matter density. The ratio F,/F, = 
o? < 1 is consistent with the discussion about the in- 
medium pion velocity. One can also show that the correc- 
tions at O(p°) do not spoil the validity of the Gell-Mann-— 
Oakes—Renner relation, see also (400), in particular both 
F,, and the quark condensate decrease with increasing 
density. For a more detailed discussion of 2-, 3- and 4- 
point functions in the medium, we reder to Ref. (404). 

The missing ingedrient in these calculations are the ef- 
fects of multi-nucleon interactions. It has recently been 
shown how these can be included in the path integral for- 
mulation (412; 413). For that, one introduces heavy fields 
H that couple to nucleon bilinears in appropriate spin- 
isospin combinations, L[U;v,a,s,p] > L[U;H;,a,s, pl, 
and letting the mass of the H-fields tend to infinity, 
cf. Fig. 3. In that way, one can formally integrate out 
the multi-fermion interactions from the generating func- 
tional. 

Further progress has been made by Girlanda et al. 
in Ref. (414). They developed a generalization of in- 
medium CHPT for finite systems. This provides a frame- 
work to study pion-nuclear bound states, for which the 
finite volume and the surface of the nucleus are impor- 
tant ingredients. The corresponding chiral counting is 
applied to the underlying pion-nucleon interactions and 
also to the relevant nuclear matrix elements. The cen- 
tral object of this approach are the Greens function in 
the presence of a nucleus, G4(X — Y) that describe the 
general process A+ X — A+Y, where X,Y represent 
some number of external pions and photons and A is a 
nucleus made of A nucleons. The presence of the nucleus 
is parameterized in terms of proton and neutron distri- 
bution functions that are taken from phenomenology. In 
the limit of uniform desnity, this approach reduce to the 
in-medium CHPT described above. As an example, the 
pion-nucleus optical potential is calculated at NLO, 


UE. D =[asnt-O [O(E7,0.2) +00], 
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where £ is the pion energy, q,@ are the out-going and 
the in-comimg three-momenta, in order, and p, is the 
proton charge density. The ellipsis in Eq. (3.7) stand 
for the contributions from the hard virtual photons and 
from the strong interaction, for details see (414). As par- 
ticularly stressed in (415), this approach allows one to 


U(E:7,¢,2) = (3.7) 
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FIG. 37 Upper panel: one- and two pion exchange diagrams 
contributing to the energy per particle at two and three loops. 
Solid and dashed lines denote nucleons and pions, respec- 
tively. Lower panel: Three-body diagrams related to 27- 
exchange with single delta excitations (double lines). These 
represent interactions between three nucleons in the Fermi 
sea. 


identify unambiguously the nuclear finite size effects and 
to disentangle the S-, P - and D-wave contributions to 
the optical potential without invoking the local density 
approximation. For a more detailed discussion concern- 
ing also the comparison with more traditional approaches 
to pion-nuclues physics, we refer to Ref. (414). 


B. Perturbative chiral nuclear dynamics 


A somewhat different path, that has turned out to 
be very successfull phenomenologically, has been taken 
by the Munich group (416-418) (for a related work, see 
Ref. (419)). Its key element is the separation of long- and 
short-distance dynamics. The ordering scheme counts 
pion massses and momentum, the Fermi momentum and 
the nucleon-delta mass splitting as small quantities 7°, 
motivated by the fact that at nuclear matter density 
kp ~ 2M,. Therefore, pions must be included and prop- 
agation effects of the delta will be resolved. The long- 
distance physics is calculated perturbatively including 
one- and two pion exchange Hartree and Fock graphs, 
see Fig.37 for some typical diagrams contributing to the 
energy per particle in nuclear matter. The short-distance 
dynamics is also treated perturbatively, either by fine- 
tuning of an UV cutoff or adjusting the parameters that 
appear at a given observable at a given order. This dif- 
fers from the treatment of the contact interactions in free 
NN scattering Within this scheme, one can reproduce the 
empirical saturation point of nuclear matter by adjusting 
(fine-tuning) one parameter related to the short-distance 
dynamics. This corresponds to a novel mechanism for nu- 
clear matter saturation due to the repulsive contribution 
to the energy per particle generated by Pauli-blocking in 
second order (iterated) one-pion exchange. This can be 
better understood by considering the realistic parame- 


#8 We unify here the two formulations with and without explicit 
deltas presented by this group. 


terization of the energy per particle in isospin symmetric 
nuclear matter, 


E(kp) =—— -a 


10m m2 (38) 


where a, are dimensionless parameters. If a,@ < 0, 
a saturation minimum will be obtained. In the deltaless 
theory, one obtains exactly such a form in the chiral limit 
and the parameters can be given in closed form. In fact, 
the experssion for ( is parameter-free. Fine-tuning the 
short-distance contribution to a, one obtains the proper 
binding energy of nuclear matter. Qualitatively, this pic- 
ture is not changed when the effects related to the finite 
pion mass and the delta-excitation are included. Further- 
more, one obtains as a by-product a realistic value of the 
nuclear matter compressibility. 

This approach has been applied and extended in 
varoius ways. Spin-orbit interactions in nuclei and hyper- 
nuclei were considered in (420-423). This led to a nice 
explanation of the very strong spin-orbit interactions in 
ordinary nuclei contrasted to the remarkably weak spin- 
orbit splitting in A hyper-nuclei. Corrections to the in- 
medium chiral condensate beyond the linear density ap- 
proximation were calculated in Refs. (424; 425). Fur- 
ther, a systematic analysis of the nuclear energy density 
functional based on a unfication of chiral pion nuclear 
dynamics with strong scalar and vector mean fields was 
performed and applied to the the properties of nuclear 
matter and finite nuclei, see Refs. (426-429). For more 
details on these interesting calculations, the reader is re- 
ferred to the original articles. 


C. EFT for halo nuclei 


A special class of nuclear systems exhibiting univer- 
sal behavior are halo nuclei. A halo nucleus consists of 
a tightly bound core surrounded by one or more loosely 
bound valence nucleons. The valence nucleons are char- 
acterized by a very low separation energy compared to 
those in the core. As a consequence, the radius of the 
halo nucleus is large compared to the radius of the core. 
A trivial example is the deuteron, which can be consid- 
ered a 2-body halo nucleus. The root mean square ra- 
dius of the deuteron (r2)'/? ~ 3 fm is about four times 
larger than the size of the constituent nucleons. Halo nu- 
clei with two valence nucleons are particularly interest- 
ing examples of 3-body systems. If none of the 2-body 
subsystems are bound, they are called Borromean halo 
nuclei. This name is derived from the heraldic symbol 
of the Borromeo family of Italy, which consists of three 
rings interlocked in such way that if any one of the rings 
is removed the other two separate. The most carefully 
studied Borromean halo nuclei are He and !!Li, which 
have two weakly bound valence neutrons (430; 431). In 
the case of °He, the core is a He nucleus, which is also 
known as the a particle. The two-neutron separation 
energy for °He is about 1 MeV, small compared to the 
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binding energy of the a particle which is about 28 MeV. 
The neutron-a (na) system has no bound states and the 
°He nucleus is therefore Borromean. There is, however, 
a strong P-wave resonance in the J = 3/2 channel of na 
scattering which is sometimes referred to as °He. This 
resonance is responsible for the binding of °He. Thus °He 
can be interpreted as a bound state of an a-particle and 
two neutrons, both of which are in P3/2 configurations. 


Because of the separation of scales in halo nuclei, they 
can be described by extensions of the pionless EFT. One 
can assume the core to be structureless and treats the 
nucleus as a few-body system of the core and the valence 
nucleons. Corrections from the structure of the core ap- 
pear in higher orders and can be included in perturbation 
theory. Cluster models of halo nuclei then appear as lead- 
ing order approximations in the EFT. A new facet is the 
appearance of resonances as in the neutron-alpha system 
which leads to a more complicated singularity structure 
and renormalization compared to the few-nucleon system 
discussed above (432). 


The first application of effective field theory methods 
to halo nuclei was carried out in Refs. (432; 433), where 
the na system (“°He”) was considered. It was found 
that for resonant P-wave interactions both the scatter- 
ing length and effective range have to be resummed at 
leading order. At threshold, however, only one combi- 
nation of coupling constants is fine-tuned and the EFT 
becomes perturbative. More recent studies have focused 
on the consistent inclusion of the Coulomb interation in 
two-body halo nuclei such as the pa and aa systems 
(434; 435). In particular, the aa system shows a sur- 
prising amount of fine-tuning between the strong and 
electromagnetic interaction. It can be understood in an 
expansion around the limit where, when electromagnetic 
interactions are turned off, the °Be ground state is ex- 
actly at threshold and exhibits conformal invariance. In 
this scenario, the Hoyle state in '1*C would appear as a 
remnant of an excited Efimov state (115). In order to bet- 
ter understand the modification of the Efimov spectrum 
and limit cycles by long-range interactions such as the 
Coulomb interaction, a one dimensional inverse square 
potential supplemented with a Coulomb interaction was 
investigated in (436). The results indicate that the coun- 
terterm required to renormalize the inverse square po- 
tential alone is sufficient to renormalize the full prob- 
lem. However, the breaking of the discrete scale invari- 
ance through the Coulomb interaction leads to a modi- 
fied bound state spectrum. The shallow bound states are 
strongly influenced by the Coulomb interaction while the 
deep bound states are dominated by the inverse square 
potential. 

Three-body halo nuclei composed of a core and two va- 
lence neutrons are of particular interest due to the possi- 
bility of these systems to display the Efimov effect (83). 
Since the scattering length can not easily be varied in 
halo nuclei, one has to look for excited states. Such stud- 
ies have previously been carried out in cluster models and 
the renormalized zero-range model (437-439). A com- 
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FIG. 38 The various one- and two-body matter density form 
factors with leading order error bands for the ground state 
of *°C in the low-energy region: Fnn(k”) black solid line; 
Fne(k?) black dotted line; Fn (k?) lighter (red) dashed line; 
F.(k?) lighter (red) dot-dashed line. 


prehensive study of S-wave halo nuclei in EFT including 
structure calculations with error estimates was recently 
carried out in Ref. (440). Currently, the only possible 
candidate for an excited Efimov state is 29°C, which con- 
sists of a core nucleus with spin and parity quantum num- 
bers J? = 0+ and two valence neutrons. The nucleus !°C 
is expected to have a 4 state near threshold, implying 
a shallow neutron-core bound state and therefore a large 
neutron-core scattering length. The value of the °C en- 
ergy, however, is not known well enough to make a defi- 
nite statement about the appearance of an excited state 
in 29°C. The matter form factors of halo nuclei can also 
be calculated in the halo EFT. As an example, we show 
the various one- and two-body matter density form fac- 
tors with leading order error bands for the ground state 
of ?°C at low momentum transfers in Fig. 38: Fnn(k?), 
Frne(k?), Fn(k?), and F.(k?). A definition of the form 
factors can be found in (440). The theory breaks down 
for momentum transfers of the order of the pion-mass 
squared (k? ~ 0.5 fm~?) where the one-pion exchange 
interaction cannot be approximated by short-range con- 
tact interactions anymore. From the slope of the matter 
form factors one can extract the corresponding radii: 


F() 1-28 lr?) 40... 


- (3.9) 


Information on these radii has been extracted from exper- 
iment for some halo nuclei. For the neutron-neutron ra- 
dius of the Borromean halo nucleus !Be for example, the 
leading order halo EFT result is \/(r2,,) = 4.1 40.5 fm. 
The value \/(r2,,)eap = 5.41.0 fm was obtained from 
three-body correlations in the dissociation of ‘Be using 
a technique based on intensity interferometry and Dalitz 
plots (441). Within the errors there is good agreement 
between both values. However, one should also keep in 
mind that there is some model dependence in the exper- 
imental result. Results for further halo nuclei are given 
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in Ref. (440). A few recent studies have also investigated 
scattering observables. In particular, in Refs. (442; 443) 
the trajectory of the possible 2°C excited state was ex- 
tended into the scattering region in order to find a reso- 
nance in n-!°C scattering. 

The simplest strange halo nucleus is the hypertriton, a 
3-body bound state of a proton, neutron, and a strange 
particle called the A. The total binding energy is only 
about 2.4 MeV. The hypertriton is not Borromean, be- 
cause the proton-neutron subsystem has a bound state, 
the deuteron. The separation energy for the A, E, = 
0.134 0.05 MeV, is tiny compared to the binding energy 
Bp = 2.224 MeV of the deuteron. The hypertriton can 
therefore also be considered a 2-body halo nucleus. It 
has been studied in both 2-body and 3-body approaches 
(444-446). A study of the hypertriton in the halo EFT 
was carried out in Ref. (364). An important feature of the 
halo EFT is the possibility to quantify theoretical errors 
through error bands. Calculations can be improved sys- 
tematically through the inclusion of higher order terms. 

Another interesting application of this effective theory 
will be the study of Coulomb excitation data from exist- 
ing and future facilities with exotic beams (such as FAIR 
and FRIB). In these experiments a nuclear beam scatters 
off the Coulomb field of a heavy nucleus. Such processes 
can populate excited states of the projectile which sub- 
sequently decay, leading to its “Coulomb dissociation” 
(447). Effective theories offer a systematic framework for 
a full quantum-mechanical treatment of these reactions. 
In summary, with new improved experimental data for 
these weakly bound nuclei, much more knowledge can 
be obtained about the structure of these interesting sys- 
tems as well as discovering whether they show universal 
behavior and excited Efimov states. 


D. Viow k potentials: construction and applications 


Nuclear interactions, like all interactions, depend on 
the resolution scale. For a momentum cutoff A, only 
details of the interactions at distances larger than 1/A 
can be resolved. The interaction potential V consists of 
2-, 3-, and higher-body terms and can be written as 

V (A) = Va(A) + V3(A) + Va(A) +... - (3.10) 
While V depends on A, observables are independent of 
A. This property can be used to construct so-called 
low-momentum potentials with a cutoff A’ < A that 
describe low-energy physics in terms of low-energy de- 
grees of freedom only. Various methods for construct- 
ing such low-momentum potentials are known. (See, e.g. 
Ref. (448) for a review.) The first construction of a low- 
momentum potential was carried out in (449; 450) based 
on the Okubo method (451; 452). Bogner and collab- 
orators have pushed this idea further and constructed 
low-momentum potentials for various realistic nucleon- 
nucleon interactions using renormalization group (RG) 
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FIG. 39 The diagonal matrix elements of the nucleon-nucleon potential in the 'So partial wave. Left panel: Various realistic 
nucleon-nucleon potentials and the chiral N?LO potential by EM (271). Right panel: Same potentials evolved down to a cutoff 


A = 2.2 fm~'. (Figure courtesy of A. Schwenk.) 


techniques. They showed that these potentials all con- 
verge to the same universal Viow x if the cutoff is lowered 
to low enough values (A S$ 2 fm~'). This is illustrated 
in Fig. 39 for the 'Sq partial wave. The left panel shows 
various realistic nucleon-nucleon potentials and the chiral 
N®LO potential by EM (271). The right panel shows the 
same potentials evolved down to a cutoff A = 2.2 fm7!. 
At this cutoff all potentials have collapsed to the same 
universal curve. 


As the resolution scale is lowered the physics previously 
present in high-momentum modes now appears in many- 
body forces that are generated through the RG trans- 
formation. The low-momentum potential constructed 
this way is phase-equivalent by construction. In most 
early calculations, the many-body forces generated by 
the RG have been neglected for simplicity. In this case 
phase equivalence only holds in the 2-body subspace. Re- 
cent advances based on similarity RG techniques, how- 
ever, suggest that these limitations can be overcome soon 
(453; 454). Low-momentum potentials have also been 
constructed for hyperon-nucleon interactions (455; 456). 
Here the various realistic potentials are less constrained 
by data and the Viow x interactions show only convergence 
in some channels. 


In Ref. (124), the three- and four-nucleon systems were 
studied using the Vow, potential supplemented by the 
leading order chiral 3N forces. This procedure was moti- 
vated by the expectation that the many body-forces have 
to reduce to the leading order chiral forces at low enough 
momentum. The free parameters in the chiral 3N force 
were then fitted to experiment for each value of the cut- 
off. If the chiral 3N force is left out, the cutoff variation 
generates the Tjon line. Since the RG evolution leaves 
the two-body observables unchanged by construction, the 
variation can only go along the Tjon line. This is in con- 
junction with the findings in the pionless EFT (120). For 
a study of potential problems with using inconsistent 2N 
and 3N low-momentum potentials in scattering calcula- 


tions see Ref. (457). 

Because of the low cutoff, the low-momentum potential 
has advantages in nuclear structure calculations where 
smaller model spaces are desirable because of the com- 
putational effort involved. For asummary of recent appli- 
cations, see Ref. (458). As the RG evolution shifts con- 
tributions between the potential and the integrals over 
intermediate states in loop integrals which are restricted 
by A, the RG transformation can eliminate sources of 
non-perturbative behavior such as strong short-range re- 
pulsion or tensor-forces (459; 460). This suggests that 
perturbative nuclear matter calculations are possible. At 
these low resolution scales, nuclear matter saturation 
would be largely driven by three-body forces. Moreover, 
perturbative nuclear matter calculations would also pro- 
vide a solid basis for the construction of a universal den- 
sity functional for nuclei with controlled errors (461; 462). 


E. Lattice simulations of many-nucleon systems 


Nuclear matter studies utilizing lattice simulations 
were pioneered by Brockmann and Frank (463), who 
calculated the quantum corrections to the Walecka 
model (464) and by Miiller et al. (465), who investi- 
gated nuclear matter properties utilizing a Hamiltonian 
that accomodates on-site and next-to-neighbour parts of 
the central, spin- and isospin-exchange nucleon-nucleon 
interactions. The first connection between chiral EFT 
and the properties of nucleon and neutron matter using 
Monte-Carlo methods was done in the groundbreaking 
work of Borasoy, Lee and Schafer (379). They laid out 
the framework for nuclear lattice simulations with chiral 
EFT and presented leading order results for hot neutron 
matter at temperatures T = 20 — 40K and densities be- 
low twice the nuclear matter density. Neutron matter 
in a periodic box based on the lattice representation of 
the chiral NLO potential (see sect. II.G) was performed 
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FIG. 40 Results for Eo /Efee versus the Fermi momentum 
kr at LO (open triangles) and NLO (filled triangles). For 
comparison, we also display the results of FP 1981 (469), 
APR 1998 (470), CMPR v6 and v8’ 2003 (471), SP 2005 (472) 
and GC 2007 (473). 


in Ref. (466), probing the density range from 2% to 8% 
of normal nuclear matter density. Dilute neutron mat- 
ter is a particularly good testing ground for chiral EFT 
applied to many-nucleon systems because of the Pauli 
suppression of three-body forces. Furthermore, neutron 
matter at kp ~ 80MeV, with kp = (3n?N)'/3/L (for 
N neutrons in a box of volume L°), is close to the so- 
called unitary limit. In this limit, the scattering length 
is infinite and the range of the interaction is zero, so 
that the scattering amplitude takes its largest possible 
value (as given by unitarity). In this limit, the only di- 
mensionful parameter describing the ground state of the 
many-fermion system is the particle density. Thus, the 
ground state energy Eo of the system obeys the simple 
relation 

Eo=tES* , (3.11) 
where € is a dimensionless measurable constant and Ef 
the ground state energy of a free Fermi gas. Due to 
its universal nature, the unitary limit can be studied 
in ultracold atomic systems like °Li or *°Ka utilizing 
Feshbach-resonance techniques. Recently measured val- 
ues for € scatter considerably and have sizeable error bars 
(for a review, see (467)). There also have been numer- 
ous calculations of € employing very different many-body 
techniques, see Ref. (398) for a recent review. 

Recent EFT simulations at LO and NLO with an im- 
proved action for 8, 12 and 16 neutrons boxes of length 
L = 10,12 and 14fm are shown in Fig. 40 (468). The 
chiral EFT results are consistent with most earlier cal- 
culations based on different methods (note the markedly 
different slope obtained in (472)). A good fit to the lat- 
tice data is obtained by (the structure of the correction 
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terms is discussed in detail in Ref. (466)) 


EK 
oo aa 40.16 kero — (0.51fm?)k% , (3.12) 
0 


with € ~ 0.31 and & ~ 0.81. This is consistent with 
the Monte Carlo studies of many-fermion systems in 
Ref. (474) but smaller than the value for € obtained 
in fixed-node Greens function Monte Carlo calculations, 
see e.g. (471). This suggests that the upper bound on 
the ground state energy in that type of approach might 
be lowered further by a more optimal fermionic nodal 
surface. Clearly, such numerical simulations of many- 
nucleon systems have become a valuable tool to fur- 
ther constrain the nuclear equation of state at moderate 
densities and lead to further insight into the physics of 
strongly-coupled many-body systems. In particular, they 
provide another nice link between nuclear and atomic 
physics as already discussed in sec. II.A. 


IV. SUMMARY AND PERSPECTIVES 


In this review, we have described the theory that has 
emerged by applying effective field theory methods to the 
nuclear force problem. This method allows for a system- 
atic derivation of nuclear forces with a direct connection 
to QCD via its symmetries. The review focused on the 
derivation of the forces and their application in the few- 
nucleon problem where most work has been carried out so 
far. However, there are many frontiers where future work 
is required. These include a better understanding of non- 
perturbative renormalization and improved renormaliza- 
tion schemes, the consistent inclusion of electroweak cur- 
rents, and the development of a consistent EFT for the 
nuclear many-body problem. The application of new 
techniques and advanced calculational methods for the 
many-body problem will be decisive to achieve the latter. 
Promising approaches include the renormalization group, 
nuclear lattice calculations, coupled cluster approaches, 
the no-core shell model, and density functional theory. 
For very low-energy processes, these approaches can be 
complemented by the pionless or halo EFT which is an 
ideal tool to unravel universal properties and establish 
connections to other fields of physics. 
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